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PROGRAM  INFORMATION 


Description  of  Program 

CBNTBM,  called  X0065  in  the  Conversationally  Oriented  Real-Time  Program- 
Generating  System  (CORPS)  Library,  is  a  computer  program  for  the  analysis  of 
plane  frame  structures  which  are  an  assemblage  of  two-dimensional  beam-column 
elements  with  linear  and/or  nonlinear  spring  supports.  The  centroidal  axis 
of  the  structure  is  assumed  to  be  a  continuous  tortuous  line  in  a  single 
plane.  CBNTBM  uses  the  finite  element  method  of  analysis  to  analyze  the 
frame-spring  system.  CBNTBM  is  intended  to  be  an  easy-to-use  program  incor¬ 
porating  the  capabilities  required  by  a  diverse  group  of  users. 

Coding  and  Data  Format 

CBNTBM  is  written  in  FORTRAN  and  is  operational  on  the  following  systems: 
a.  WES  Honeywell  DPS/1. 

_b.  Harris  500  computers  which  are  located  at  most  district  Corps  offices. 

c^.  Control  Data  Corporation,  Cybernet  Computer  Service's  CDC  CYBER 
systems. 

Data  can  be  input  either  interactively  at  execute  time  or  from  a  prepared 
data  file  with  line  numbers.  Output  may  be  directed  to  an  output  file  or 
come  directly  back  to  the  terminal. 

How  to  Use  CBNTBM 


A  short  description  of  how  to  access  the  program  on  each  of  the  three  systems 
is  provided  below.  It  is  assumed  that  the  user  knows  how  to  sign  on  the 
appropriate  system  before  trying  to  use  CBNTBM.  In  the  example  initiation  of 
execution  commands  below,  all  user  responses  are  underlined,  and  each  should 
be  followed  by  a  carriage  return. 

WES  Honeywell  System 

After  the  user  has  signed  on  the  system,  the  two  system  commands  FORT  and  NEW 
get  the  user  to  the  level  to  execute  the  program.  Next,  the  user  issues  the 
run  command 

RUN  WESLIB/CORPS/X0065.R 


to  initiate  execution  of  the  program.  The  program  is  then  run  as  described 
in  this  user* 8  guide.  The  data  file  should  be  prepared  prior  to  issuing  the 
RUN  command.  An  example  initiation  of  execution  is  as  follows,  assuming  a 
data  file  had  previously  been  prepared: 


HIS  SERIES  600  ON  03/04/81  AT  13.301  CHANNEL  5647 

USER  ID  -  ROKACASECON 

PASSWORD  -  XXXXXXXXXXXXXX 

SYSTEM?  FORT  NEW 

READY 

*RUN  WESLIB/CORPS/X0065 , R 


CDC,  Cybernet  System 


The  log-on  procedure  is  followed  by  a  call  to  the  CORPS  procedure  file 
OLD,  CORPS/UN=CECELB 

to  access  the  CORPS  Library.  The  file  name  of  the  program  is  used  in  the 
command 

BEGIN,, CORPS, X006 5 

to  initiate  execution  of  the  program.  An  example  is: 

S4/01/25  10.32.51  AC2E5DA 

EASTERN  CYBERNET  CENTER  SN487  NOS  1 . 4/531 . 281/20AD 

FAMILY:  KOE 

USER  NAME:  CER0C2 

PASSWORD  - 

XXXXXXX 

TERMINAL:  510.NAMIAF 

RECOVER/CHARGE:  CHARGE,  CEXXXXX,  YYYYYY 
/OLD,  CQRPS/UN=CECELB 
/BEGIN,, CORPS/X0065 


Local  District  Harris  Systems 


After  the  user  has  signed  on  the  system,  the  command  to  execute  the  CORPS 
program  will  be 

*CORPS,X0065 

An  example  to  illustrate  the  log-on  and  execution  procedure  on  one  Harris  500 
is  shown  below.  There  may  be  some  differences  at  some  local  Corps  sites. 

"ACOE  -  VICKSBURG" 

USER  #?  NNNNWES  WESXXX 

**  Good  Morning  25  Jan  84  9:56:31 

VED  HARRIS  500 


*CORPS ,X0065 


b 


How  to  Use  CORPS 


The  CORPS  system  contains  many  other  useful  programs  which  may  be  catalogued 
from  CORPS  by  use  of  the  LIST  command.  The  execute  command  for  CORPS  on  the 
WES  system  is: 

RUN  WESLIB/ CORPS/CORPS , R 

ENTER  COMMAND  (HELP, LIST, BRIEF, MESSAGE, EXECUTE,  OR  STOP) 

*?LIST 


On  the  Cybernet  system,  the  commands  are: 

OLD,CORPS/UN=CECELB 
BEGIN , , CORPS , CORPS 

ENTER  COMMAND  (HELP, LIST, BRIEF, MESSAGE, EXECUTE,  OR  STOP) 
*?LIST 


On  the  Harris  local  systems,  the  commands  are: 
*CORPS 

ARE  YOU  USING  A  PRINTER  TERMINAL  OR  CRT? 

ENTER  P  OR  C 
C 

ENTER  COMMAND  (BRIEF, EXECUTE, LIST, HELP, STOP) 
LIST 
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This  user's  guide  documents  a  computer  program  called  CBNTBM  that  can  be 
used  for  analysis  of  a  plane  frame  structure  which  is  an  assemblage  of  two- 
dimensional  beam-column  elements  with  linear  and/or  nonlinear  spring  supports. 
The  work  in  writing  the  computer  program  and  the  user's  guide  was  accomplished 
with  funds  provided  to  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  Vicksburg,  Miss.,  under  the  Civil  Works  Research  and  Development  Pro¬ 
gram  of  the  Office,  Chief  of  Engineers,  U.  S.  Army  (OCE),  as  part  of  the  Struc¬ 
tural  Engineering  Research  Program  work  unit  of  the  Soil-Structure  Interaction 
(SSI)  Studies  Project. 

The  computer  program  and  user's  guide  were  written  by  Dr.  William  P. 
Dawkins,  P.E.,  of  Stillwater,  Okla . ,  under  Contract  No.  DACW39-81-M-0715  with 
WES. 

Dr.  N.  Radhakrishnan,  Special  Technical  Assistant,  Automatic  Data 
Processing  (ADP)  Center,  WES,  and  SSI  Studies  Project  Manager,  coordinated 
and  monitored  the  work.  Messrs.  H.  Wayne  Jones  and  Reed  L.  Mosher,  Computer- 
Aided  Design  Group,  provided  technical  assistance  in  developing  and  evaluating 
the  program.  Mr.  Donald  R.  Dressier,  Structures  Division,  Civil  Works 
Directorate,  was  the  OCE  point  of  contact. 

Commanders  and  Directors  of  WES  during  the  development  of  the  program 
were  COL  Nelson  P.  Conover,  CE,  and  COL  Tilford  C.  Creel,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  NON- SI  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (metric) 
units  as  follows: 


Multipl 


To  Obtain 


feet 

0.3048 

meters 

inches 

0.0254 

meters 

kip  (force)-feet 

0.00135582 

newton-meters 

kips  (1000  lb  force) 

4.448222 

kilonewtons 

kips  (force)  per  foot 

0.0145939 

newtons  per  meter 

kips  (force)  per  inch 

0.17512685 

newtons  per  meter 

kips  (force)  per  square  inch 

6.894757 

m<  iscals 

pound  (force) -inches 

0.1129848 

i  ton-meters 

pounds  (force) 

4.448222 

tons 

pounds  (force)  per  foot 

14.5939 

w  is  per  meter 

pounds  (force)  per  inch 

175.12685 

newtons  per  meter 

pounds  (force)  per  square  inch 

6.894757 

kilopascals 

square  inches 

0.00064516 

square  meters 

■»  wv> 
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USER'S  GUIDE:  COMPUTER  PROGRAM  FOR  ANALYSIS  OF 
TWO-DIMENSIONAL  BEAM-C01UMNS  WITH  NONLINEAR 
SUPPORTS  (CBNTBM) 


PART  1:  INTRODUCTION 
Background 

1.  A  beam-column  is  usually  understood  to  be  an  initially  straight 
structural  member  in  which  axial  internal  forces  amplify  the  effects  produced 
by  lateral  loads.  A  finite  element  model  and  a  general-purpose  computer  pro¬ 
gram  for  analysis  of  straight  beam-columns--CBEAMC--are  described  in  Dawkins 
(1982) .* ** 

2.  This  report  presents  an  extension  of  the  beam-column  analysis  to  a 
special  class  of  plane  frame  structure  in  which  the  frame  is  composed  of  an 
assemblage  of  beam-column  elements.  Documentation  for  a  computer  program-- 
CBNTBM-- employing  the  technique  is  provided."*5' 

Organization  of  Report 


3.  The  remainder  of  this  report  is  organized  as  follows: 

a.  Part  II:  Describes  the  structural  system  considered  and  the 

mathematical  model  used  for  analysis. 

b.  Part  III:  Presents  the  force-displacement  relationships  for  the 
mathematical  model  and  describes  the  computational  procedure  used 
for  solution. 

c.  Part  IV:  Describes  the  computer  program. 

d.  Part  V:  Guide  for  data  input. 

e.  Part  VI:  Provides  example  solutions  obtained  with  the  program. 


*  W.  P.  Dawkins.  1982.  "User's  Guide:  Computer  Program  for  Analysis  of 
Beam-Column  Structures  with  Nonlinear  Supports  (CBEAMC),"  Instruction  Re¬ 
port  K-82-6,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg, 
Miss . 

**  CBNTBM  is  designated  X0065  in  the  Conversationally  Oriented  Real-Time 

Program-Generating  System  (CORPS)  library.  Three  sheets  entitled  "PROGRAM 
INFORMATION"  have  been  hand-inserted  inside  the  front  cover  of  this  report. 
They  present  general  information  on  CBNTBM  and  describe  how  it  can  be 
accessed.  If  procedures  used  to  access  this  and  other  CORPS  programs  should 
change,  recipients  of  this  report  will  be  furnished  revised  versions  of  the 
"PROGRAM  INFORMATION"  sheets. 


Remarks 


4.  CBNTBM  has  been  checked  for  computational  accuracy  within  reasonable 
limits.  However,  there  may  exist  unusual  situations  which  were  not  anticipated 
which  may  cause  the  program  to  produce  questionable  results.  It  is  the  respon¬ 
sibility  of  the  user  to  use  good  engineering  judgment  to  determine  the  validity 
of  the  results. 


PART  II:  STRUCTURAL  SYSTEM  ANT)  FINITE  ELEMENT  MODEL 


5.  The  structure  is  assumed  to  be  a  special  class  of  plane  frame  struc¬ 
ture  in  which  the  centroidal  axis  is  a  continuous  tortuous  line  in  a  single 
plane.  Example  tortuous  line  structures  are  shown  in  Figure  1.  The  procedure 
presented  in  this  report  is  restricted  to  those  structures  in  which  the 


f 


e.  Arch  f.  U- frame 

Figure  1.  Example  tortuous  line  structures 


eentroidal  axis  is  composed  of  a  sequence  of  initially  straight  segments.  Ap¬ 
proximate  solutions  for  initially  curved  structures  may  be  obtained  by  replac¬ 
ing  the  curved  structure  with  straight  segments  which  are  chords  of  the  origi¬ 
nal  curve. 


Global  Coordinate  System  and  Nodes 


6.  Several  coordinate  systems  are  required  to  describe  the  geometry  and 
to  interpret  the  behavior  of  the  two-dimensional  structure.  The  eentroidal 
axis  of  the  structure  is  assumed  to  lie  in  the  x-y  plane  of  a  right-handed 
Cartesian  x-y-z  coordinate  system  referred  to  as  the  global  coordinate  system. 
The  origin  and  orientation  of  the  global  system  are  arbitrary. 

7.  The  geometry  of  the  structure  is  established  by  the  global  x-y 
coordinates  of  nodes  (joints)  on  the  eentroidal  axis.  Nodes  are  described  as 
either  "dictated"  or  "interne  ate." 


Dictated  Nodes 

8.  Dictated  nodes  are  those  required  to  establish  the  overall  geometry 
of  the  structure  at: 

a.  Each  end  of  the  structure. 

b.  The  juncture  of  two  straight  segments  of  the  eentroidal  axis. 

c.  Each  change  in  cross-sectional  or  material  properties  ^--mod¬ 
ulus  of  elasticity;  I — moment  of  inertia;  A--cross-sectional 
area) . 

d.  The  point  of  application  of  a  concentrated  load. 

e.  The  point  of  application  of  a  fixed  support  or  concentrated 

spring  (linear  or  nonlinear)  support. 

f.  The  beginning  and  end  of  a  distributed  load. 

g.  The  beginning  and  end  of  distributed  spring  (either  linear 

or  nonlinear)  supports. 


Intermediate  Nodes 

9.  Intermediate  nodes  are  those  defined  between  dictated  nodes  as  nec¬ 
essary  to  improve  accuracy  of  the  solution  (discussed  later)  or  to  provide  in 
creased  detail  of  the  variation  of  structural  response. 


Node  Numbers 


10.  Nodes  are  identified  by  integer  numbers  starting  with  node  1  at 
one  end  of  the  structure  and  proceeding  sequentially  along  the  centroidal 
axis.  Only  the  node  numbers  and  global  coordinates  of  the  dictated  nodes  are 
required  as  input  to  the  computer  program.  Global  coordinates  (and  node 
numbers)  for  intermediate  nodes  are  generated  at  equal  intervals  (along  a 
straight  line  or  circular  arc;  see  Part  V,  "Guide  for  Data  Input")  between 
successive  pairs  of  dictated  nodes. 

Element  Identification 

11.  An  element  is  defined  as  the  straight  segment  of  the  structure  be¬ 
tween  successive  nodes  and  is  identified  by  the  lower  of  its  two  terminal  node 
numbers;  i.e.,  element  i  lies  between  nodes  i  and  i+1.  The  structure  is  thus 
composed  of  n  nodes  (node  1  at  one  end,  node  n  at  the  other)  and  n-1  elements. 

Local  Coordinate  Systems 

12.  Local  coordinate  systems  are  established  for  convenience  in  describ 
ing  behavior  or  interpreting  the  response  of  loads,  spring  supports,  internal 
forces,  and  displacements.  The  local  coordinate  system  associated  with  each 
item  is  described  below. 

Displacements 

13.  Displacements  of  the  system  are  completely  defined  by  global  compo¬ 
nents  of  the  nodal  translations  u  (parallel 
to  the  global  x  axis)  and  v  (parallel  to 
the  global  y  axis),  and  rotation  6  about 
an  axis  perpendicular  to  the  x-y  plane  (see 
Figure  2). 

14.  All  displacements  are  assumed  to 
be  "small";  i.e.,  they  do  not  alter  the  basic 
geometry  of  the  structure.  Rotations  0  are 
small  enough  to  permit  the  approximation 


y (Global ) 

1  v 


u 

i 

—  x  (Global ) 


cos  6=1  and  sin  6=6. 


Figure  2.  Nodal 
displacements 


y  (Global ) 


x  (Global ) 


Figure  3.  Applied 
concentrated  loads 


Applied  Concentrated  Loads 


15.  Concentrated  loads  may  be  applied  to  any 
node.  These  loads  are  assumed  to  be  applied  as  com¬ 
ponents  parallel  to  the  global  x  and  y  axes  (F  and 
Fy  ,  respectively)  or  as  couples  C  about  axes 
perpendicular  to  the  x-y  plane  (see  Figure  3). 


Fixed  Supports 


16.  A  fixed  support  is  an  external  influence  which  results  in  a  spe¬ 
cific  value  (either  zero  or  nonzero)  of  one  or  more  of  the  global  components 
(u  ,  v  ,  0)  of  displacement  at  a  node  regardless  of  any  other  effects. 


Concentrated  Linear  Springs 


17.  Concentrated  linear  springs  produce  forces  on  the  structure  which 
are  directly  proportional  to  the  displacements  of  the  node  to  which  the 
springs  are  attached: 

a.  A  translational  spring  and  its  attendant  local  coordinate 
system  are  shown  in  Figure  4.  The  local  axis  extends  from 


Figure  4.  Linear  concentrated  spring 


the  attachment  node  toward  the  fixed  base  of  the  spring  making 
an  angle  a  (positive  counterclockwise)  with  the  global  x 
axis.  The  deformation  of  the  spring  may  be  expressed  in  terms 
of  the  nodal  translations  as 

A  =  u.  cos  a  +  v.  sin  a 
s  x  1  (1 


The  force  produced  on  the  node  is 


where  k  is  the  spring  stiffness  and  F  acts  on  the  node 
s  s 

parallel  to  xg  and  opposite  in  sense  to  A  . 

A  rotational  spring  (Figure  5)  resists  rotation  of  the  node. 

y  (Global ) 


Figure  5.  Linear  rotation 
spring 
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The  moment  produced  on  the  node  by  the  rotational  spring  is 


C  =  R  0. 
s  si 


where  R  is  the  spring  stiffness  and  C  acts  opposite  to  the 
s  s 

nodal  rotation  0.  . 

l 


Concentrated  Nonlinear  Translational  Springs 


18.  A  concentrated  nonlinear  translational  spring  (Figure  6)  produces  a 
force  which  resists  translation  of  the  node  to  which  it  is  attached.  The  re¬ 
sisting  force  is  a  function  of,  but  not  directly  proportional  to,  the  trans¬ 
lation  of  the  node.  The  resistance-deformation  function  of  the  spring  shown 
schematically  in  Figure  6b  must  be  provided  as  input  to  CBNTBM  (Part  V) . 
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Node  i 


x  (Global) 


Nonlinear  concentrated  spring  b.  Force-deformation  relationship 

Figure  6.  Nonlinear  spring  characteristics 


The  local  coordinate  system  associated  with  the  nonlinear  spring  is  shown  in 
Figure  6a.  The  deformation  of  the  spring  is 


A  =  u.  cos  a  +  v.  sin  a  (4) 

s  1  1 

The  resisting  force  produced  by  the  spring  is  obtained  from  the 

resistance-deformation  function.  F  acts  parallel  to  the  local  x  axis 

s  s 

and  opposite  in  sense  to  the  deformation  A  .  (Note:  Nonlinear  rotational 

s 

springs  are  not  accommodated  by  the  CBNTBM.) 

Element  Coordinate  System 

19.  Element  force-displacement  relationships  as  well  as  the  effects  of 
all  distributed  loads  and  distributed  springs  (either  linear  or  nonlinear)  are 
established  in  a  local  coordinate  system  peculiar  to  each  element.  A  local 
coordinate  system  and  typical  element  are  shown  in  Figures  7  and  8.  As  stated 
previously,  the  element  is  identified  by  the  lower  of  its  two  terminal  node 
numbers.  The  local  xi  axis  is  positive  from  node  i  to  node  i+1.  The  local  z. 
axis  is  parallel  to  (and  in  the  same  positive  sense  as)  the  global  z  axis. 

The  local  y^  axis  is  established  so  that  the  x^-y^-z^  system  forms  a  right- 
handed  Cartesian  system. 

20.  The  following  assumptions  are  employed  in  interpreting  input  data 
provided  to  CBNTBM: 

a.  Distributed  loads  are  described  as  extending  between  input  pairs 
of  dictated  nodes.  If  the  terminal  nodes  of  any  element  (e.g., 
element  i  between  nodes  i  and  i+1)  are  within  the  specified 
distribution,  the  distributed  load  is  interpreted  as  being 
parallel  to  the  local  axes  (either  x^  or  y..)  for  element  i,  re¬ 
gardless  of  whether  the  terminal  nodes  of  the  distribution  and 
the  terminal  nodes  of  element  i  fall  along  a  single  straight 
line. 

b.  Distributed  springs  are  described  as  extending  between  input 

pairs  of  dictated  nodes.  If  the  terminal  nodes  of  any  element 
are  within  the  specified  distribution,  the  distributed  spring 
is  interpreted  as  resisting  displacement  components  parallel  to 
the  local  axes  (x^  or  y^)  for  element  i,  regardless  of  whether 

the  terminal  nodes  of  the  distribution  and  the  terminal  nodes 
of  element  i  fall  along  a  single  straight  line. 
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Figure  8.  Typical  element  in  local  coordinate  system 


21.  The  force-deformation  relationship  for  a  nonlinear  spring  (either 
concentrated  or  distributed)  is  assumed  to  be  provided  as  a  piecewise  linear 
curve,  as  shown  in  Figure  9.  For  any  deformaton  A  parallel  to  the  spring 


Figure  9.  Characteristics  of  nonlinear  springs 


line  of  action,  the  total  resisting  force  (concentrated  or  distributed)  may 
be  expressed  as 

Q  (or  q)  =  Qq  (or  qQ)  +  kA  (5 

Thus,  although  an  iterative  solution  is  required,  during  any  iteration  a  non¬ 
linear  spring  may  be  replaced  by  a  combination  of  a  fixed  load  (Qq  or  qQ)  and 
a  linear  spring  with  stiffness  k  . 


PART  III:  FORCE-DISPLACEMENT  RELATIONSHIPS  AND  SOLUTION  PROCEDURE 


xA 


8* 

8 
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Element  Force-Displacement  Relationships 

22.  An  element  is  defined  as  that  portion  of  the  structure  between 
adjacent  nodes.  For  reference  in  the  following  paragraphs,  an  element  is  iden¬ 
tified  by  the  node  number  at  its  left  end  as  shown  in  Figure  7. 

23.  Each  element  is  characterized  by  modulus  of  elasticity  E  ,  cross- 
sectional  area  A  ,  moment  of  inertia  I  ,  and  length  h  and  is  subjected  to 
distributed  loads  and  springs  as  shown  in  Figure  8. 

24.  Interaction  of  internal  axial  stress  resultants  with  lateral  de¬ 
flections  will  affect  the  bending  response  of  the  element.  When  this  effect  is 
included,  the  bending  resistance  of  the  element  is  developed  using  a  constant 
average  axial  stress  resultant  given  by 

p  =  r  <ui*i '  ui>  <6) 

25.  The  relationship  between  element  end  forces  and  end  nodal  displace¬ 
ments,  obtained  by  procedures  given  in  standard  texts,*  is 

=  +  -  <7) 

where 

f1  =  If1  f1  1T  =  [f1  .  f1  .  M1  f1  f1  M1  1  T 

~  |~i  ~i+lj  I  x,i  y,i  i  x,i+l  y,i+l  l+lj 

=6xi  vector  of  element  end  forces 

Sg  =  6  x  6  axial  and  beam  bending  stiffness  matrix  including  the  effects 
of  the  axial  stress  resultant  on  bending 

=  6  x  6  stiffness  matrix  representing  the  effects  of  distributed 
springs 

u1  =  Tu.  u.  1T  =  fu.  v.  e.  u.A1  v...  e.^l1 

~  l~i  ~i+l|  I  l  l  l  l+l  l+l  i+l| 

=6xi  vector  of  end  nodal  displacements 

*  R.  W.  Clough  and  J.  Penzien.  1975.  Dynamics  of  Structures,  McGraw-Hill, 

New  York. 

J.  S.  Przmieniecki .  1968.  Theory  of  Matrix  Structural  Analysis,  McGraw- 

Hill,  New  York. 

0.  C.  Zienkiewicz.  1971.  The  Finite  Element  Method  in  Engineering  Science, 
McGraw-Hill,  New  York. 
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If1  .  f1  ...1T  =  ff1  .  f1  .  M1  . 

l_~e ,  1  ~e ,  l+lj  Lex.x  ey.i  e,x 

f1  f1  M1  1 
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=  6  x  1  vector  of  fixed  end  forces  due  to  distributed  loads 


26.  Coefficients  for  matrices  S*  ,  S*  ,  and  f*  are  given  in  Fig¬ 


ures  10-12  and  Table  1. 


Solution  Procedure 


27.  After  all  local  force-displacement  relations  (elements  and  concen¬ 
trated  springs)  are  transformed  to  the  global  coordinate  system,  satisfaction 
of  joint  equilibrium  results  in  a  system  of  3n  (for  a  structure  with  n  nodes) 
simultaneous  equations  which  may  be  solved  for  the  nodal  displacements.  Back 
substitution  of  nodal  displacements  into  force-displacement  relations  allows 
evaluation  of  all  internal  forces. 


Iterations 


28.  Whenever  nonlinear  spring  supports  are  present  and/or  the  effect  of 
the  axial  stress  resultant  on  bending  is  included,  an  iterative  solution  of 
the  simultaneous  equations  is  required.  The  iterative  solution  is  initiated 
by  evaluating  axial  stress  resultants  and  nonlinear  spring  stiffnesses  for 
zero  displacements.  On  each  succeeding  iteration,  the  solution  of  the  preced¬ 
ing  iteration  is  used  for  evaluating  system  properties.  The  process  is  con¬ 
tinued  until  the  results  of  two  successive  iterations  differ  by  an  acceptably 
small  amount. 


Effect  of  Node  Spacing  on  Solution 


29.  The  finite  element  procedure  described  above  produces  "exact"  solu¬ 
tions  for  systems  possessing  the  following  characteristics: 

a.  The  structure  is  piecewise  prismatic  and  composed  only  of 
straight  segments. 

Only  concentrated  external  loads  and  linearly  varying  lateral 
(local  y)  loads  are  present.  Or,  if  the  effect  of  the  axial 
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Figure  12.  Element  fixed  end  forces  due  to  distributed  loads 


stress  resultant  on  element  bending  stiffness  is  excluded, 
distributed  (local  x)  loads  vary  linearly. 

c.  Only  concentrated  springs  are  present. 

d.  Adjacent  elements  are  approximately  equal  in  length.  Adjacent 
elements  having  substantially  different  lengths  may  result  in 
significant  round-off  errors  in  the  solution  of  simultaneous 
equations . 

30.  The  number  of  nodes  (and  elements)  used  in  the  finite  element  model 
may  affect  the  accuracy  of  the  solution.  In  general,  as  the  number  of  nodes 
(and  elements)  is  increased,  the  solution  tends  to  converge  to  the  "exact" 
solution.  There  is  no  rule  of  thumb  to  provide  the  necessary  number  of  nodes 
for  acceptable  results.  It  may  be  necessary  to  perform  several  solutions  for 
various  numbers  of  nodes  and  node  spacings  to  ensure  that  an  adequate  solution 
has  been  obtained. 


PART  IV:  COMPUTER  PROGRAM 


31.  A  computer  program  implementing  the  analytical  process  described 
above  has  been  written  in  FORTRAN  language  for  interactive  execution.  With 
minor  revisions,  the  program  will  be  operational  on  any  computer  employing  a 
60-bit  (  15-decimal-digit)  word  length.  For  systems  using  fewer  than  15  deci¬ 
mal  digits,  it  may  be  necessary  to  perform  some  arithmetic  operations  in 
double-precision. 

Input  Data 

32.  Input  data  may  be  provided  from  the  user  terminal  during  execution 
or  from  a  predefined  data  file.  CBNTBM  has  been  written  to  generate  automat¬ 
ically  intermediate  data  values  from  a  minimum  of  user  input  data.  Details 
of  required  user  input  data  are  presented  in  Part  V.  The  processes  and 
assumptions  used  in  converting  user  input  data  to  intermediate  data  values  are 
described  below. 

Displacement  and  Load  Sign  Conventions 


33.  Positive  directions  for  displacements  and  loads  are  as  follows: 

a.  Nodal  displacements  u  and  v  are  positive  if  the  node 
translates  in  the  positive  global  x  and  y  directions, 
respectively. 

b.  Nodal  rotation  is  positive  if  the  node  rotates  counterclockwise. 

c.  Concentrated  nodal  loads  F x  and  F  are  positive  if  they  act 

in  the  positive  global  x  and  y  directions,  respectively.  Dis¬ 
tributed  loads  and  q^  are  positive  if  they  act  in  the 

positive  element  (local)  x  and  y  directions,  respectively. 


Distributed  Nonlinear  Springs 

34.  Distributed  nonlinear  springs  produce  distributed  forces  propor¬ 
tional  to  the  translation  displacements.  The  deformation  of  the  spring  is 
equal  to  and  has  the  same  sign  as  the  displacement.  The  distributed  force 
produced  on  an  element  by  a  distributed  nonlinear  spring  is  positive  if  the 
force  acts  in  the  positive  element  local  x  and  y  coordinate  directions. 


Distributed  spring  characteristics  are  assumed  to  be  described  by  pairs  of  re¬ 
sisting  force-deformation  coordinate  values  similar  to  those  for  concentrated 
springs . 

Data  Generation 

35.  Concentrated  applied  loads,  concentrated  springs,  and  fixed  sup¬ 
ports  occur  at  the  nodes  on  the  structure.  All  other  data  (cross-section 
properties  E  ,  A  ,  and  I  and  distributed  loads  and  springs)  are  distrib¬ 
uted  over  a  range  of  nodes.  The  user  is  required  to  provide  the  locations 
(node  numbers)  at  which  concentrated  effects  occur  and  the  terminal  values 
(beginning  and  ending  node  numbers)  of  distributed  quantities.  These  define 
the  "dictated"  nodes  in  the  model.  The  program  also  provides  for  defining 
nodes  intermediate  to  the  dictated  nodes.  The  following  paragraphs  describe 
the  procedures  and  assumptions  used  to  obtain  global  coordinates  and  system 
characteristics  at  the  intermediate  nodes. 

Node  Generation 

36.  The  user  describes  the  basic  geometry  of  the  structure  by  providing 
node  numbers  and  global  coordinates  of  the  "dictated"  nodes  beginning  with 
node  1,  proceeding  in  increasing  order  of  node  number,  and  ending  with  the 
data  for  the  last  node  (node  n)  in  the  structure.  Nodes  with  numbers  between 
consecutive  pairs  of  "dictated"  nodes  are  the  "intermediate"  nodes.  The  user 
also  designates  whether  the  segment  of  the  structure  between  consecutive 
dictated  nodes  is  a  straight  line  or  a  circular  arc.  Coordinates  of  the  inter¬ 
mediate  nodes  are  generated  automatically  as  follows.  (The  number  of  elements 
in  a  segment  is  equal  to  the  difference  between  the  consecutive  pairs  of  input 
dictated  node  numbers . ) 

a.  A  straight  segment  is  divided  into  elements  of  equal  length 
with  the  length  of  each  element  equal  to  the  length  of  the 
segment  divided  by  the  number  of  elements  in  the  segment. 

b.  If  a  circular  arc  segment  is  specified,  the  global  coordinates 
of  the  center  of  the  arc  are  also  required  as  input.  Global 
coordinates  of  the  intermediate  nodes  are  generated  around  the 
circular  arc  at  equal  spacing.  The  curved  segment  is  approxi¬ 
mated  by  a  polygonal  arc  composed  of  straight  elements  between 
nodes.  The  accuracy  of  this  approximation  depends  on  the  number 
of  nodes  in  the  arc. 


37.  Node  data  may  be  provided  at  any  location  on  the  structure  even 
though  such  data  are  not  associated  with  a  "dictated"  node.  These  "voluntary 
dictated"  nodes  may  be  used  to  control  the  spacing  of  intermediate  nodes;  i.e 
element  length. 

38.  Data  for  some  "dictated"  nodes  may  be  omitted  if  the  automatic  gen 
eration  of  node  data  described  above  will  result  in  a  node  at  the  required 
location;  e.g.,  at  the  beginning  and/or  end  of  a  distributed  load  or  spring 
or  at  the  location  of  a  concentrated  load,  concentrated  spring,  or  fixed 
support . 

Cross-Section  Properties 

39.  Cross-section  properties  are  provided  at  the  terminals  of  a  range 
of  node  numbers.  The  modulus  of  elasticity  E  must  be  constant  within  the 
range.  Cross-section  area  A  and  moment  of  inertia  I  are  assumed  to  vary 
linearly  between  values  provided  at  the  terminals.  This  linear  variation  is 
then  approximated  by  a  sequence  of  prismatic  elements  with  the  properties  of 
each  element  determined  as  shown  in  Figure  13.  Nonzero  values  of  E  ,  A  , 
and  I  must  be  provided  (or  generated)  for  every  element  in  the  structure. 

If  an  element  is  encompassed  by  more  than  one  range  of  input  values,  the 
properties  calculated  for  each  range  are  cumulative. 

Distributed  Loads  and  Springs 

40.  Distributed  loads  and  springs  are  described  by  values  at  the  termi¬ 
nals  of  a  range  of  node  numbers.  Distributed  loads  and  springs  are  implicitly 
associated  with  the  local  coordinate  system  for  the  element  on  which  they  act. 
For  example,  if  a  distributed  load  q^  is  specified  over  the  range  of  nodes  j 
to  m  in  Figure  13,  the  load  will  be  interpreted  as  acting  perpendicular  to 
each  element  in  the  range.  A  linear  variation  of  effects  (loads  or  springs) 
between  values  at  the  terminals  is  assumed,  and  values  at  nodes  within  the 
range  are  determined  as  shown  in  Figure  13. 

Output  Data 

41.  Output  data  may  be  directed  to  the  user  terminal,  to  a  data  file, 
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Figure  13.  Linear  interpolation  of  input  data 


or  to  both.  The  user  selects  the  type  and  extent  of  output  from  the  sections 
described  below. 


Echoprint  of  Input  Data 

42.  This  section  contains  a  tabulation  of  all  input  data  currently 
available  to  the  program.  Output  of  this  section  is  optional. 


Node  Coordinates  and  Displacements 

43.  A  tabulation  of  the  global  coordinates  and  global  components  of  the 
displacements  for  all  nodes  in  the  structure  is  available.  The  user  may  spec¬ 
ify  a  complete  tabulation,  or  select  one  or  more  ranges  of  nodes  for  output, 
or  omit  output  of  this  section. 

Element  End  Forces 

44.  A  tabulation  of  the  forces,  in  the  element  coordinate  system,  on  the 
ends  of  each  element  in  the  structure  is  available.  The  user  may  specify  tne 
extent  of  output  for  this  section. 

Reactions  at  Fixed  Supports 

45.  A  tabulation  of  the  global  components  of  reactions  at  all  fixed 
supports  is  available. 

Forces  in  Concentrated  Linear  Springs 

46.  A  tabulation  of  spring  deformation  and  force  in  the  local  coordi¬ 
nate  system  is  available  for  all  concentrated  linear  springs. 

Forces  in  Distributed  Linear  Springs 

47.  A  tabulation  of  the  cumulative  force  in  the  distributed  linear 
springs  at  each  end  of  each  element  in  the  element  coordinate  system  is  avail¬ 
able.  The  user  may  specify  the  extent  of  output  for  this  section. 

Forces  in  Concentrated  Nonlinear  Springs 


48.  A  tabulation  of  the  deformation  and  force  in  each  nonlinear  concen¬ 
trated  spring  is  available. 

Forces  in  Distributed  Nonlinear  Springs 


49.  A  tabulation  of  the  cumulative  force  in  the  distributed  nonlinear 


springs  at  each  end  of  each  element  in  the  element  coordinate  system  is  avail 
able.  The  user  may  specify  the  extent  of  output  for  this  section. 


Output  Sign  Conventions 


50.  All  forces  and  translational  displacements  are  positive  if  they 
act  in  the  associated  positive  coordinate  direction  (either  global  or  local) 
All  counterclockwide  rotations  and  moments  are  positive. 


PART  V:  GUIDE  FOR  DATA  INPUT 


Source  of  Input 

51.  Input  data  may  be  supplied  from  a  predefined  data  file  or  from  the 
user  terminal  during  execution.  If  data  are  supplied  from  the  user  terminal, 
prompting  messages  are  printed  to  indicate  the  amount  and  character  of  data 
to  be  entered. 


Data  Editing 

52.  When  all  data  for  a  problem  have  been  entered,  the  user  is  offered 
the  opportunity  to  review  an  echoprint  of  the  currently  available  input  data 
and  to  revise  any  or  all  sections  of  the  input  data  before  execution  is  at¬ 
tempted.  When  data  are  edited  during  execution,  each  section  must  be  entered 
in  its  entirety. 


Input  Data  File  Generation 


53.  After  data  have  been  entered  from  the  terminal,  either  initially  or 
after  editing,  the  user  may  direct  the  program  to  write  the  input  data  to  a 
permanent  file  in  input  data  file  format. 


Data  Format 


54.  All  input  data  (whether  supplied  from  the  user  terminal  or  from  a 
file)  are  read  in  free-field  format: 

a.  Data  items  must  be  separated  by  one  or  more  blanks  (comma 
separators  are  not  permitted). 

b.  Integer  numbers  must  be  of  the  form  NNNN 

c.  Real  numbers  may  be  of  the  form 

txxxx,  ixx.xx,  or  ±xx.xxE+ee. 

d.  User  responses  to  all  requests  for  control  by  the  program 
for  alphanumeric  input  may  be  abbreviated  by  the  first  let¬ 
ter  of  the  indicated  work  response;  e.g., 

ENTER  'YES'  OR  'NO' --respond  Y  or  N 

ENTER  'CONTINUE'  or  'END' —  respond  C  or  E 


Carriage  return  responses  alone  will  cause  abnormal  termina¬ 
tion  of  the  program. 


Sections  of  Input 

55.  Input  data  are  divided  into  the  following  sections: 

a.  Heading. 

b.  Geometry  Data. 

c.  Cross-Section  Properties. 

d.  Applied  Loads. 

e.  Fixed  Supports. 

f.  Linear  Spring  Data. 

g.  Nonlinear  Spring  Data. 

h.  Termination. 

56.  When  data  are  entered  from  the  terminal,  data  sections  a  through  g 
may  be  input  in  any  order.  When  data  are  entered  from  a  predefined  data  file, 
the  Heading  section  (a)  must  be  entered  first;  other  sections  (b  through  g) 
may  be  entered  in  any  order. 

57.  When  data  are  entered  from  the  terminal,  the  user  is  cued  for 
Termination.  Termination  for  a  data  file  is  discussed  later. 

Minimum  Required  Data 

58.  Data  sections  a,  b,  and  c  are  always  required.  At  least  one  of  sec 
tions  e,  f,  and  g  is  required.  It  is  the  responsibility  of  the  user  to  ensure 
that  sufficient  supports  (either  fixed  or  spring  supports)  are  provided  to  in¬ 
hibit  all  rigid  body  motions  of  the  system;  i.e.,  to  ensure  a  stable  structure 

Units 

59.  The  program  recognizes  the  following  units: 

Inches  Feet  Pounds  Kips 

Default  is  to  inches  and  pounds.  Angular  units  are  either  degrees  or  radians 
as  explained  below. 

60.  Each  data  section  may  be  entered  with  any  combination  of  units 
for  length  and  force  as  specified  by  the  user. 


Predefined  Data  File 


61.  In  addition  to  the  general  format  requirements  given  in  para¬ 
graph  4,  above,  the  following  pertain  to  a  predefined  data  file  and  to  the 
input  data  description  which  follows: 

a.  Each  line  must  commence  with  a  nonzero,  positive  line  number, 
denoted  LN  below. 

b.  A  line  of  input  may  require  both  alphanumeric  and  numeric  data 
items.  Alphanumeric  data  items  are  enclosed  in  single  quotes 
in  the  following  paragraphs. 

c.  A  line  of  input  may  require  a  keyword.  The  acceptable  abbrevi¬ 
ation  for  the  keyword  is  indicated  by  underlined  capital  letters 
e.g.,  the  acceptable  abbreviation  for  the  keywork  ’ Linear '  is 
'Ll*. 

d.  Lower  case  words  in  single  quotes  indicate  a  choice  of  keywords 
defined  following. 

e.  Items  designated  by  uppercase  letters  and  numbers  without  quotes 
indicate  numeric  data  values.  Numeric  data  values  are  either 
real  or  integer  according  to  standard  FORTRAN  variable  naming 
conventions . 

f.  Data  items  enclosed  in  brackets  [  ]  may  not  be  required.  Data 
items  enclosed  in  braces  {  }  indicate  special  note  follows. 

g.  Input  data  are  divided  into  the  sections  discussed  in  para¬ 
graph  5  above.  Except  for  the  Heading,  each  section  consists 
of  a  header  line  and  one  or  more  data  lines.  The  header  line 
serves  the  multiple  purposes  of:  indicating  the  end  of  the 
preceding  section;  identifying  the  data  section  to  follow;  and 
indicating  the  units  associated  with  the  section. 

h.  Comment  lines  may  be  inserted  in  the  input  file  by  enclosing 
the  line,  following  the  line  number,  in  parentheses.  Comment 
lines  are  ignored;  e.g., 

1234|6  (THIS  LINE  IS  IGNORED). 

Input  Description 

62.  Heading — One  (1)  to  four  (4)  lines 

a.  Contents 

LN  'heading' 

b.  Definition 

'heading'  =  any  alphanumeric  information  (70  characters 
maximum) 

c.  Discussion--a  header  line  must  begin  with  a  single  quote  if 
the  line  begins  with  the  first  letter  of  any  of  the  section 
titles  or  keywords  described  below. 


63.  Geometry  Data--Three  (3)  to  twenty-two  (22)  lines 

a.  Control — One  (1)  line 
a.l  Contents 

LN  'Geometry'  [’units'] 

a.  2  Definitions 

'Geometry'  =  section  title 

'Units'  =  'Inches'  or  'Feet' 

'Inches'  assumed  if  omitted 

b.  Node  Coordinate  Data--Two  (2)  to  twenty-one  (21)  lines 

b. l  Contents 

LN  NODE  X  Y  [’segment  type'  XC  YC] 
b.2  Definitions 

NODE  =  node  number 
X,Y  =  global  coordinates 

['segment  type']  =  'Straight'  or  'Circular';  'Straight' 

assumed  if  omitted 

[XC,  YC]  =  global  coordinates  of  center  of  'Cir¬ 
cular'  arc  segment 

b.3  Discussion 

b.3.a  Data  for  NODE  =  1  must  be  provided  first. 
b.3.b  Node  numbers  on  successive  lines  are  assumed  to 
indicate  the  ends  of  segments  of  the  structure. 
b.3.c  Missing  nodes  are  generated  at  equal  intervals 

between  input  node  numbers  along  a  straight  line 
or  along  a  circular  arc  according  to  the  'segment 
type'  data  provided. 

b.3.d  The  last  node  provided  (denoted  NLAST  below) 
determines  the  end  of  the  structure. 

64.  Cross-Section  Properties  Data--Two  (2)  to  twenty-two  (22)  lines 
a.  Control — One  (1)  line 

a.l  Contents 

LN  'Properties'  ['units'] 
a. 2  Definitions 

'Properties’  =  section  title 

['units']  =  'Inches'  or  'Feet'  and  'Pounds'  or  'Kips' 
'Inches'  and  'Pounds'  assumed  if  omitted 
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b.  Properties  Data--One  (1)  to  twenty-one  (21)  lines 
b.l  Contents 

LN  N0DE1  NODE 2  E  A1  SI1  [A2  SI 2] 
b.2  Definitions 

NODE1,  N0DE2  =  node  numbers  at  beginning  and  end  of 
range,  respectively 

E  =  modulus  of  elasticity  between  N0DE1  and 
NODE  2 

A1  =  cross-section  area  at  N0DE1 
SI1  =  cross-section  moment  of  inertia  at  N0DE1 
(A2,  SI2]  =  cross-section  area  and  moment  of  inertia, 
respectively,  at  N0DE2;  assumed  equal  to 
values  at  N0DE1  if  omitted 

b.3  Discussion 

b.3.a  Properties  for  elements  between  N0DE1  and  N0DE2 
are  determined  by  linear  interpolation  from  sec¬ 
tion  properties  specified  by  the  end  values. 
b.3.b  N0DE1  and  N0DE2  must  satisfy: 

1<  N0DE1  <  NODE 2  <  NLAST 

b.3.c  Properties  for  elements  specified  on  more  than  one 
line  are  cumulative. 

65.  Applied  Loads — Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines;  en¬ 
tire  section  may  be  omitted 

a.  Control — One  (1)  line 
a.l  Contents 

[LN  ' LOads ’  ['units]] 

a.  2  Definitions 

'LOads'  =  section  title 

['units']  =  'Inches'  or  'Feet'  and  'Pounds'  or  'Kips'; 

'Inches'  and  'Pounds'  assumed  if  omitted 

b.  Data  Lines  for  Corcentrated  Loads 

b. l  Contents 

[LN  'Concentrated'  NODE  FX  FY  C] 
b.2  Definitions 

’Concentrated'  =  keyword 

NODE  =  node  number  at  concentrated  load 


FX,  FY  =  global  x  and  y  components  of  concen¬ 
trated  load  at  NODE 

C  =  applied  couple  at  NODE  (positive  CCW) 

c.  Data  Lines  for  Distributed  Loads 
c.l  Contents 

[LN  'Distributed'  'direction'  N0DE1  Q1  N0DE2  [Q2]] 

c. 2  Definitions 

'Distributed'  =  keyword 

'direction'  =  'X'  if  distributed  load  acts  in  LOCAL 

x  coordinate  direction  for  elements  be¬ 
tween  N1  and  N2 

=  'Y'  if  load  acts  in  LOCAL  y  coordinate 
direction 

N0DE1,  N0DE2  =  node  numbers  at  start  and  end,  respec¬ 
tively,  of  distribution 
Q1  =  magnitude  of  distributed  load  at  N0DE1 
[Q2]  =  magnitude  of  distributed  load  at  NODE  2; 
assumed  equal  to  Q1  if  omitted 

d.  Discussion 

d. l  Data  may  not  be  described  for  nodes/elements  beyond  the 

limits  of  the  structure;  i.e., 

for  concentrated  loads:  1  <  NODE  <  NLAST 
for  distributed  loads:  1  <  N0DE1  <  N0DE2  <  NLAST 

d.2  Magnitudes  of  distributed  loads  at  nodes  between  N0DE1 
and  N0DE2  are  determined  by  linear  interpolation  from 
values  at  the  ends  of  the  distribution. 
d.3  Values  for  nodes/elements  specified  on  more  than  one  line 
are  cumulative. 

66.  Fixed  Supports — Zero  (0)  or  two  (2)  to  eleven  (11)  lines;  entire 
section  may  be  omitted 

a.  Control — One  (1)  line 
a.l  Contents 

[LN  'Fixed'  ['units']] 
a. 2  Definitions 

'Fixed'  =  keyword 


CKS 


V  V  V  % 
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['units']  =  'Inches'  or  'Feet';  'Inches'  assumed  if 
omitted 

b.  Data  Lines — One  (1)  to  ten  (10)  lines 
b.l  Contents 

[LN  NODE  XD  YD  R] 
b.2  Definitions 

NODE  =  node  number  at  specified  conditions 

XD  =  displacement  in  X  GLOBAL  coordinate  direction  to 
be  enforced  at  NODE 

=  'Free'  if  node  is  unrestrained  in  X  direction 
YD  =  displacement  in  Y  GLOBAL  coordinate  direction 
=  'Free'  if  node  is  unrestrained  in  Y  direction 
R  =  rotation  in  radians  to  be  enforced  at  NODE  (posi 
tive  CCW) 

=  'Free'  if  NODE  is  free  to  rotate 
b.3  Discussion 

b.3.a  Data  may  not  be  described  for  nodes  beyond  the 
limits  of  the  structure;  i.e., 

1  <  NODE  <  NLAST 

b.3.b  Only  one  set  of  specified  displacements  may  be 
imposed  at  a  single  node. 

67.  Linear  Spring  Data--Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines 
entire  section  may  be  omitted 

a.  Control — One  (1)  line 
a.l  Contents 

[LN  ' Linear'  ['units']] 

a.  2  Definitions 

'Linear'  =  section  title 

[’units']  =  'Inches'  or  'Feet'  and  'Pounds'  or  'Kips'; 

'Inches'  and  'Pounds'  assumed  if  omitted 

b.  Data  Lines  for  Concentrated  Linear  Springs 

b. l  Contents 

[LN  'Concentrated'  NODE  ANGLE  S  R] 


b.2  Definitions 

'Concentrated'  =  keyword 


NODE  =  node  number  at  spring 
ANGLE  =  angle  in  degrees  between  GLOBAL  X  axis 
and  line  of  action  of  translation  re¬ 
sisting  spring  (see  Figure  4) 

S  =  stiffness  (force/unit-deformation)  of 
translation  spring 

R  =  stiffness  (moment/radian)  of  rotation 
resisting  spring 

Data  Lines  for  Distributed  Linear  Springs 
c.l  Contents 

[LN  'Distributed'  'direction'  N0DE1  SI  N0DE2  [S2] ] 

c. 2  Definitions 

'Distributed'  =  keyword 

'direction'  =  'X'  if  spring  resists  displacements  in 
local  (element)  x  direction 
=  'Y'  if  spring  resists  displacements  in 
local  (element)  y  direction 
N0DE1,  N0DE2  =  node  numbers  at  start  and  end,  respec¬ 
tively,  of  distribution 
SI  =  spring  stiffness  (force/unit-length/ 
unit-displacement)  at  start  of 
distribution 

[S2]  =  spring  stiffness  at  end  of  distribution; 
assumed  equal  to  SI  if  omitted 

Discussion 

d. l  Data  may  not  be  described  for  nodes/elements  beyond  the 

limits  of  the  structure;  i.e., 

for  concentrated  springs:  1  <  NODE  <  NLAST 
for  distributed  springs:  1  <  N0DE1  <  N0DE2  <  NLAST 

d.2  Distributed  springs  resist  displacements  in  the  local 
(element)  coordinate  directi  is. 
d.3  Distributed  spring  stiffnesses  at  nodes  between  N0DE1  and 
N0DE2  are  determined  by  linear  interpolation  from  values 
at  the  ends  of  the  distribution. 
d.4  Values  for  nodes  specified  on  more  than  one  line  are 
cumulative . 


Nonlinear  Spring  Data--Zero  (0)  or  four  (4)  to  sixty-four  (64) 
entire  section  may  be  omitted 

a.  Header — One  (1)  line 
a.l  Contents 

ILN  'NONlinear'  [’units']] 

a.  2  Definitions 

'NONlinear'  =  section  title 

['units']  =  'Inches'  or  'Feet'  and  'Pounds'  or  'Kips’; 

'Inches'  and  'Pounds'  assumed  if  omitted 

b.  Data  Lines  for  Concentrated  Nonlinear  Springs--Three  (3) 
lines  for  each  concentrated  spring 

b. l  Line  1  contents 

[LN  'Concentrated'  NODE  ANGLE  SPTS  DMUL  FMUL] 

Line  2  contents 

[LN  DEF(l)  DEF(2)  .  .  .  DEF(NPTS)] 

Line  3  contents 

[LN  FORCE ( 1 )  FORCE (2)  .  .  .  FORCE (NPTS)] 
b.2  Definitions 

'Concentrated'  =  keyword 

NODE  =  node  number  at  spring 
ANGLE  =  angle  (in  degrees)  between  GLOBAL  x  axis 
and  spring  line  of  action  (see  Figure  5) 
NPTS  =  number  of  points  on  resisting  force- 
deformation  curve;  minimum  of  two  (2) 
points  required;  maximum  of  eight  (8) 
points  permitted 

DMUL  =  scale  factor  for  curve  deformation 

coordination  values;  must  be  positive 
and  nonzero 

FMUL  =  scale  factor  for  curve  resisting  force 
coordinate  values 

DEF(  )  =  curve  deformation  coordinate  value;  NPTS 
values  on  iine;  must  proceed 
sequentially  from  most  negative  to  most 
positive 


FORCE (  )  =  curve  resisting  force  coordinate  value 
corresponding  to  DEF(  ) 

Data  Lines  for  Distributed  Nonlinear  Springs — Minimum  of  Two 
(2)  groups  of  three  (3)  lines  for  each  distribution 
c.l.a  Group  1,  Line  1,  Contents 

(LN  'Distributed'  'direction'  N0DE1  NPTS  DMUL 
FMUL] 

Group  1,  Line  2,  Contents 

[LN  DEF(l)  DEF(2)  .  .  .  DEF(NPTS) ] 

Group  1,  Line  3,  Contents 
[LN  FORCE ( 1 )  FORCE (2)  .  .  .  FORCE (NPTS)] 
c.l.b  Definitions  for  Group  1 

'Distributed'  =  keyword  indicating  beginning  of 
distribution 

'direction*  =  'X'  or  'Y'  indicating  distributed 

spring  resists  x  or  y  displacements  in 
local  (element)  coordinate  directions 
N0DE1  =  node  at  beginning  of  distribution 
NPTS  =  number  of  points  provided  on  nonlinear 
resisting  force-displacement  curve  at 
NODE  1;  minimum  of  two  (2)  points  re¬ 
quired;  maximum  of  eight  (8)  points 
permitted 

DMUL  =  scale  factor  for  curve  deformation  co¬ 
ordinate  values;  positive,  nonzero 
FMUL  =  scale  factor  for  curve  resisting  force 
coordinate  values 

DEF(  )  =  curve  deformation  coordinate  value;  NPTS 
values  on  one  line;  must  proceed  from 
most  negative  to  most  positive 
FORCE (  )  =  curve  resisting  distributed  force  cor¬ 
responding  to  DEF(  ) 
c.2.a  Group  2,  Line  1,  Contents 

[LN  'control'  NODE 2  NPTS  DMUL  FMUL] 

Group  2,  Line  2,  Contents 

[LN  DEF(l)  DEF(2)  .  .  .  DEF(NPTS) ] 


Group  2,  Line  3,  Contents 

[LN  FORCE ( 1 )  FORCE ( 2 )  .  .  .  FORCE(NPTS)] 

c. 2.b  Definitions  for  Group  2 

'control'  =  'Continue'  or  'End';  'Continue'  indicates 
distribution  continues  and  another  set  of 
Group  2  lines  follows;  'End'  indicated  this 
is  the  last  curve  in  this  distribution 
N0DE2  =  node  for  this  nonlinear  curve 
NPTS ,  DMUL ,  FMUL,  DEF(  ),  FORCE (  )  as  in  Group  1 
d.  Discussion 

d. l  Final  deformation  and  resisting  force  curve  coordinates 

are  products  DMUL*DEF(  )  and  FMUL*FORCE(  ),  respectively. 
d.2  Resisting  force-deformation  curve  must  be  single-valued 
at  every  point. 

d.3  Resisting  force-deformation  curve  is  assumed  to  vary 
linearly  between  input  points. 
d.4  Resisting  force  is  assumed  to  be  constant  at  first  or 
last  value  provided  on  curve  for  deformations  beyond 
the  range  of  deformation  coordinates  provided 
d.5  Several  concentrated  nonlinear  springs  may  be  attached 
at  a  single  node. 

d.6  Characteristics  of  distributed  nonlinear  springs  at 
intermediate  nodes  on  a  distribution  are  obtained  by 
linear  interpolation  between  adjacent  curves. 
d.7  Overlapping  distributions  result  in  cumulative  effects 
at  nodes  in  the  overlap. 

d.8  Data  may  not  be  specified  for  nodes/elements  beyond  the 
limits  of  the  structure;  i.e., 

for  concentrated  springs:  1  <  NODE  <  NLAST 
for  distributed  springs:  1  <  N0DE1  <  N0DE2  <  NLAST 

Termination — One  (1)  line 

a.  Contents 
LN  'FINish' 

b.  Definition 

'FINish'  =  keyword  indicating  end  of  input  data  for  a  problem. 
Initiates  data  checking  and  solution 


70.  Data  items  enclosed  in  brackets  [  ]  may  be  omitted.  Braces  {  } 
enclosing  data  lists  indicate  choose  one;  arrow  (-»)  indicates  default  if  item 
is  omitted. 

71.  Heading — One  (1)  to  four  (4)  lines 
LN  'heading' 

[LN  'heading'] 

[LN  'heading'] 

[LN  'heading'] 

72.  Geometry  Data--Three  (3)  to  twenty-two  (22)  lines 

a.  Control — One  (1)  line 


LN  'Geometry'  Wnches’l  h|OUD?Si 
’  'Feet  'Kips' 

b.  Data  Lines--0ne  (1)  to  twenty-one  (21)  lines 


LN  NODE  X  Y 


(-►'Straight* 
|  'Circular' 


-  'Circular'  XC  vCJJ 

73.  Cross-Section  Properties — Two  (2)  to  twenty-two  (22)  lines 

a.  Control--0ne  (1)  line 

Tu  p*’  Inches  '1  (-*'  Pounds'] 

LN  'Properties  -  - 

'Feet'  'Kips' 

_  J  ^ 

b.  Data  Lines--0ne  (1)  to  twenty-two  (22)  lines 
LN  N0DE1  &0DE2  E  A1  SI1  [A2  SI2] 

74.  Applied  Loads — Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines 


]1)  line 

Inches'! 

f-*'  Pounds  '1 

[  'Feet'  J 

\  'Kips’  •  JJ 

b.  Data  lines  for  Concentrated  Loads — One  (1)  line  for  each  loaded 
node 

[LN  'Concentrated'  NODE  FX  FY  C] 

c.  Data  Lined  for  Distributed  Loads — One  (1)  line  for  each 
distribution 

[LN  'Distributed'  N0DE1  Q1  N0DE2  [Q2] ] 

75.  Fixed  Supports — Zero  (0)  or  two  (2)  to  eleven  (11)  lines 
a.  Control — One  (1)  line 


tv  .v  Th’ Inches']  [-►'Pounds'!!' 

LN  'Fixed’  |j  ,gMt,  |  |  ’ Kips '  jj 


b.  Data  Lines--One  (1)  line  for  each  supported  node 
[LN  NODE  XD  YD  R] 

76.  Linear  Spring  Data — Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines 
a.  Control--0ne  line 


lines 


LN 

' Linear' 

(-*'  Inches '  1 

(-»'  Pounds  '1 

|  'Feet'  j 

1  ’Kips' 

b.  Data  Lines  for  Concentrated  Linear  Springs--0ne  (1)  line  for 
each  spring 

[LN  'Concentrated'  NODE  ANGLE  S  R] 

£.  Data  Lines  for  Distributed  Linear  Springs--One  (1)  line  for 
each  distribution 


[LN  'Distributed 


'  .{*} 


N0DE1  SI  N0DE2  [S2] ] 


77.  Nonlinear  Spring  Data--Zero  (0)  or  four  (4)  to  sixty-four  (64) 

a.  Control — One  (1)  line 

[» {i°r}] 

b.  Data  Lines  for  Concentrated  Nonlinear  Springs--Three  (3) 

i 

lines  for  each  spring 
Line  1 : 

[LN  'Concentrated'  NODE  ANGLE  NPTS  DMUL  FMUL] 

Line  2: 

[LN  DEF(l)  DEF ( 2 )  .  .  .  DEF(NPTS)] 

Line  3: 

[LN  FORCE ( 1 )  FORCE (2)  .  .  .  FORCE (NPTS)] 

c.  Data  Lines  for  Distributed  Nonlinear  Springs--Minimum  of  two 
(2)  groups  of  three  (3)  lines  for  each  distribution 

Group  1 ,  Line  1 : 


LN  'Distributed 


'  W 


N0DE1  NPTS  DMUL  FMUL 


Group  1,  Line  2: 

(LN  DEF(l)  DEF(2)  .  .  .  DEF (NPTS) ] 

Group  1,  Line  3: 

[LN  FORCE ( 1 )  FORCE (2)  .  .  .  FORCE (NPTS)] 


s 


T 


Figure  14.  Example  1 — gabled  frame 


1,  2,  6,  7,  and  8)  are  required  to  define  the  structural  geometry.  Addi¬ 
tional  nodes  along  the  loaded  member  were  used  to  provide  more  detail  of  the 
variation  of  internal  forces  in  this  region. 

81.  Data  were  entered  during  execution  from  the  terminal  as  shown  on 
pages  43  and  44.  An  echoprint  of  the  input  data  is  shown  on  page  45,  and  a 
listing  of  the  input  data  file  generated  by  the  program  is  given  on  page  48. 

82.  Output  data  were  printed  to  the  terminal  as  shown  on  pages  46  and 
47.  There  is  no  approximation  involved  in  this  solution,  and  the  results  are 
identical  with  those  obtained  by  other,  closed-form  solutions. 

83.  The  program  concludes  a  successful  solution  with  the  normal 
termination  message  shown  on  page  47.  Any  abnormal  termination  prior  to  this 
message  may  result  in  loss  of  data  files  generated  by  the  program. 
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Execution  of  CBNTBM  with  Data  Entered  Interactively 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE!  02/15/83  TIME:  08 ! 34 ! 43 

ARE  INPUT  DATA  TO  BE  READ  FROM  TERMINAL  OR  FILE? 

ENTER  ' TERMINAL '  OR  'FILE' 

I>T 

INFUT  DATA  SECTIONS  AND  STATUS 


SECTION 

NUMBER 

1  .  .  . 

SECTION  CONTENTS 
.......  HFADTNR  . . 

DATA 

AVAILABLE 

ACTIVE 

UNIT1 

......  .RFflMETRY  TlATA.  ........ 

. . INCHES 

3  ♦  t  * 

. ..^FTTinN  PRHPFRT  TES .  .  . 

. . INCHES 

POUND 

4  .  .  . 

» . INCHES 

POUND 

5 .  ♦  » 

. . INCHES 

6  ♦  «  * 

.  .  INCHES 

POUND' 

POUND 

. . INCHES 

ENTER  SECTION 

NUMBER  < 1  TO  7  >  FOR  DATA  TO 

BE  INPUT 

ENTER  'FINISHED'  AFTER  ALL  DATA  SECTIONS  HAVE  BEEN  INPUT 
ENTER  'STATUS'  TO  OBTAIN  ABOVE  STATUS  TABLE 

I  I 

ENTER  NUMBER  OF  HEADING  LINES  <1  TO  4) 


I  >  1 

ENTER  HEADING  LINE  1  (72  CHARACTERS  MAXIMUM) 
i:  EXAMPLE  1  -  GABLED  FRAME 

ENTER  SECTION  NUMBER  (1  TO  7)  FOR  DATA  TO  BE  INPUT 

I  >2 

GEOMETRY  DATA  --  ENTER  DATA  UNITS 
'INCHES'  OR  'FEET'  OR  'DEFAULT' 

1/  F 

ENTER  GEOMETK  '  DATA  ONE  LINE  AT  A  TIME 
ENTER  END'  IF  FINISHED  WITH  THIS  SECTION 

NODE  NODE  COORDINATES  SEGMENT  TYPE 

NO .  X  Y  ('S'  OR  ' C '  > 

(FT)  (FT) 

i:-i  0  0 
I >2  0  25 
1/  6  30  00 
1/  7  100  25 
I >8  100  0 
1/  END 


CENTER  COORDINATES 
X 

(FT)  (FT' 


(Continued) 


CO  CO  U1  tfl 


5-  Lul 


Execution  of  CBNTBM  with  Data  Entered  Interactively  (Continue 


ENTER  SECTION  NUMBER  <1  TO  7)  TOR  DATA  TO  BE  INPUT 

i:  3 

SECTION  PROPERTIES  --  ENTER  DATA  UNITS 
INCHES'  OR  'FEET'  AND  'POUNDS'  OR  'KIPS' ,  OR  DEFAULT' 

I  D 

ENTER  SECTION  PROPERTIES  ONE  LINE  AT  A  TIME 
ENTER  END'  IF  FINISHED  WITH  V HIS  SECTION 

-  - SECTION  PROPERTIES--  ----- 

NODE  NUMBER  MODULUS  OF  -> . START -  -  - STOP - 

START  STOP  ELASTICITY  AREA  INERTIA  AREA  INERTIA 

( PS I i  ( SOIN  >  (IN**4>  (SUIN'  '  I N ♦ *  4 ' 

IM  S  30.  E4  47.1  9760 
I:  END 

ENTER  SECTION  NUMBER  (t  TO  7)  FOR  DATA  TO  BE  INPUT 

I  -1 

LOAD  DATA  --  ENTER  DATA  UNITS 

'INCHES'  OR  FEET'  AND  'FOUNDS'  OR  KIPS'.  OR  'DEFAULT 

i:  f  F 

ENTER  LOAD  TYPE  --  'CONCENTRATED'  OR  DISTRIBUTED' 

ENTER  'END'  IF  FINISHED  WITH  THIS  SECTION 

i:  D 

ENTER  DISTRIBUTED  LOAD  DATA  ONE  LINE  AT  A  TIME 
ENTER  'END'  WHEN  FINISHED  WITH  DISTRIBUTED  LOADS 

LOAD  - START -  - - STOP - . 

DIRECT  NODE  LOAD  NODE  LOAD 

(  X  .  Y  }  (i  T'  ‘  P/T  • 

I  X  2  -650  6 
2  -1300  6 
ND 

ENTER  LOAD  T V F E  --  'CONCENTRATED'  OR  'DISTRIBUTED' 

ENTER  'END'  IF  FINISHED  WITH  THIS  SECTION 

I  END 

ENTER  SECTION  NUMBER  (1  TO  7)  FOR  DATA  TO  BE  INPUT 

I’j 

FIXED  SUPPORT  DATA  --  ENTER  DATA  UNITS 
INCHES'  OR  FEET  OR  'DEFAULT' 

i:d 

ENTER  FIXED  SUPPORT  DATA  ONE  LINE  AT  A  TIME 
ENTER  'HD'  IF  FINISHED  WITH  THIS  BFCTIDN 

SPECIFIED  DISPLACEMENTS 

SUPPORT  (ENTER  FREE  IF  UNSPE  IF) HO 

NODE  NO.  X-DI3F  Y-DISP  ROTATION 

•IN)  IN'-  (RAD- 

I"1  OOF 
I>8  0  0  F 
I  -END 

ENTER  SECTION  NUMBER  (1  TO  7)  FOR  DATA  TO  BE  INPUT 

I  F 

INPUT  COMPLETE 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL,  TO  A  FILE,  TO  BOTH,  OK  NEITHER7 
ENTER  'TERMINAL',  'FILE',  'BOTH  ,  OR  NEITHER' 

I  T 


(Continued) 


.  V  V  *. 


L 


i 

a 


Execution  of  CBNTBM  with  Data  Entered  Interactively  (Continued) 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE!  02/ 15/83  TIMF:  08!38:22 

I. --INPUT  DATA 

1 .--HEADING 


EXAMPLE  1  - 

GABLED  FRAME 

2. --GEOMETRY 

DATA 

NODE 

NODE  COORDINATES 

SEGMENT 

CENTER 

COORDINATES 

NO. 

X 

Y 

TYPE 

X 

Y 

(FT) 

(FT  ) 

(FT  ) 

(FT) 

1 

0.00 

0.00 

STRAIGHT 

n 

0.00 

25.00 

STRAIGHT 

6 

50.00 

50.00 

STRAIGHT 

7 

100.00 

25.00 

STRAIGHT 

8 

100.00 

0.00 

3. --SECTION  PROPERTIES 


NODE  NUMBER 
START  STOP 


4. --LOAD  DATA 


MODULUS  OF 
ELASTICITY 
<PSI  ) 
3.00E+07 


<  - SECTION  PROPERTIES - > 

<  - START-  — >  < - STOP - > 

AREA  INERTIA  AREA  INERTIA 

( SO I N  )  <IN**4)  ( SOIN )  (IN**4> 

47.10  9760,00  47.10  9760. CO 


4. A. --CONCENTRATED  LOADS 
NONE 


4. B. --DISTRIBUTED  LOADS 

LOAD  < - START - > 

DIRECT  NODE  LOAD 

NO,  <F'/FT> 

X  2  -  650 . 00 

Y  2  -1300.00 


- stop - > 

NODE  LOAD 

NO.  (P/TT) 

6  -650.00 

6  -1300.00 


--FIXED  SUPPORT  DATA 
SUPPORT 
NODE 
NO. 

1 

8 


6. --LINEAR  SPRING  DATA 
NONE 


7. --NONLINEAR  SPRING  DATA 
NONE 


SPECIFIED  DISPLACEMENTS 
X-DISP.  Y-DISP •  ROTATION 

(IN)  (IN)  (RAD) 

o.oo  o.oo  free: 

0.00  0.00  FREE 


(Continued) 


TOW 


HO  YOU  WANT  TCI  EDIT  INPUT  DATA?  ENTER  'YES'  OR  'NO' 


IN 

DO  YOU  WANT  INPUT  DATA  SAVED  IN  A  TILE?  ENTER  'YES'  OR  '  NO ' 

I  >  Y 

ENTER  FILE  NAME  FOR  SAVING  INPUT  DATA  (6  CHARACTERS  MAXIMUM) 

I !  CBNT 1 

DO  YOU  WANT  TO  CONTINUE  SOLUTION"?  ENTER  YES'  OR  'NO' 

I :  V 

DO  YOU  UANT  AXIAL  FORCE  EFFECTS  ON  PENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION-?  ENTER  'YES'  OR  'NO¬ 
TH 


SOI UTION  COMPLETE 

DO  YOU  UANT  RESUITS  WRITTEN  TO  TOUR  TERMINAL.  TO  A  FILE.  OR  BOTH? 
ENTER  'TERMINAL'.  'FILE  •  OR  '  BGTH  ' 

I  T 

FOLLOWING  OUTPUT  DATA  ARP  AVAIIAPL!' 

SECTION 


NUMBER  SECTION  CONI  ENTS 

1  . NODT  COORDINATES  AND  D  I  SRI  ACEME  NTS 

2  . ELEMENT  END  EOREF". 

3  . REACTIONS  AT  FIXED  SUE  FORTS 


COMPLETE  OUTPUT  CONTAINS  APPROXIMATELY  32  LINES. 
OUTPUT  CONTROL 

ENTER  'ALL'.  'SELECTIVE'.  FINISHED'.  OR  HE L F  ’ 

I  A 

ENTER  DESIRED  OUTPUT  UNITS. 

'INCHES'  OR  ' f  F  E  I  AND  POUNDS  OR  KIPS'.  OR  DEFAULT 

I  D 


(Continued) 


Execution  of  CBNTBM  with  Data  Entered  Interactively  (Concluded) 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2H  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
hate:  02/15/83  time:  03:40:13 

HEADING 

EXAMPLE  1  -  GABLED  FRAME 


NODE 

COORDINATES  AND 

DISPLACEMENT 

S  IN  GLOBAL 

COORDINATE 

SYSTEM 

NODE 

NODE  CO OR 

DINATES 

- NODE 

:  DISPLACEMENTS - 

NO. 

X  (IN) 

Y  (IN- 

X  (IN) 

Y  (IN) 

KOI  *  RAD  .• 

1 

0.00 

0.00 

0  . 

0  * 

4. 3856-05 

9 

0.00 

300.00 

2.S90E- 01 

-  1 . 2  9  4  E  -  0  2 

-2.673E-03 

3 

150 ♦ 00 

375 » 00 

5. 194 E -01 

-5.439E  -01 

-3.774E  03 

4 

300.00 

450.00 

7 . 4  9  6  E  -  0 1 

-1 . 0 1 2  E  +  0  0 

-  2 . 1  4 1 E  -  0  3 

5 

450.00 

525.00 

3. 106E-01 

-1 . 1 40C+00 

4.770E-04 

6 

600 ♦ 00 

600 . 00 

6.937E-01 

-9.039E-01 

2 . 3  3  4  E  -  0  3 

7 

1200.00 

300.00 

1 . 133E+00 

-4.313E-03 

-1 . 9  6  2  E - 0  3 

3 

1200.00 

0 . 00 

0. 

0. 

-  4 . 6  8  3  E  -  0  3 

EL. 

NO  . 

1 
2 

3 

4 

5 

6 

© 

*  REACTIONS  AT  FIXED  SUPPORTS  IN  GLOBAL  COORDINATE  SYSTEM 


NODE 

X-REACT 

Y-REACT 

MOM-REACT 

NO. 

(P) 

(F) 

(P-IN) 

1 

1  . 771E  +  04 

6.094E+04 

0. 

8 

-1 .771E  +  04 

2 . 031E+04 

0. 

DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS? 

ENTER  'YES'  OR  'NO'. 

1:  N 

OUTPUT  COMPLETE. 

DO  YOU  WANT  TO  EDIT  INPUT  DATA  FOR  THE  PROBLEM 
JUST  COMPLETED?  ENTER  'YES'  OR  'NO' 

I  N 

DO  YOU  UANT  TO  MAKE  ANOTHER  'CBNTBM'  RUN?  ENTER  YES'  OR  'NO' 

1:  N 

**********  NORMAL  TERMINATION  ********** 


ELEMENT  END 


X-FORCE 

(P) 

6.094E+04 
4 . 309E+04 
3.401E+04 
2.492E+04 
1 • 584E+04 
2 . 492E+04 
2 . 03 1 E+04 


FORCES  IN  ELEMENT 

- LEFT  END - 

Y-FORCE 
(P  ) 

-1.771E+04  0. 

4 . 653E+04  5. 

2 . 342E+04 
1 .025E+04 
-7.91 9E+03 
-1 .025E+04 
1 . 771E+04 


COORDINATE 


MOMENT 
<P-IN> 

312E+06 
766E+05 
219E+04 
-4.414E+06 
-1 .563E+06 
5.312E+06 


9  . 
-4  . 


-6 

-3 


SYSTEM 

- RIGHT 


X-FORCE 

<  R ) 

094E  +  04 
40 1 E  +  04 
492E+04 
-1 .534E+04 
-6 . 755E+03 
-2 . 492E+04 
-  2 . 0  3 1 E  +  0  4 


END- 


Y-FORCE 

(P> 

1 .771E+04 
•2.842E+04 
1 .025ES04 
7.919E+03 
2  •  609E  +  04 
1  >  0  2  5  E  +  0  4 
-1 .771E  +  04 


MOMENT 
<  P  - 1 N  ) 
-5.312E+06 
9 . 766E+05 
4 . 2 19E+06 
4 ■ 414C  +  06 
1 .563E+06 
-5 . 3 1 2E  +  06 
0. 


& 
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Input  Data  File  for  Example  1  Created  by  CBNTBM 


■ 

1000 

EXAMPLE  1  -  GABLED  FRAME 

1010 

GEOMETRY  F 

1020 

10.  0. 

1030 

2  0.  2 . 500E+01 

1040 

6  5.000E+01  5.000E+01 

kv 

1050 

7  1 . 000E  +  02  2  •  500E  +  0 1 

■£ 

1060 

3  1.C00E+02  0. 

1070 

PROPERTIES  I  P 

■ 

10S0 

1  S  3.000E407  47.10 

9760.00  47 

1090 

LOADS  F  F 

1100 

D  X  2  -650.00 

6  -650  *  00 

1110 

D  Y  2  -1300.00 

6  -1300.00 

»y« 

1120 

FIXED  I 

-  \ 

1130 

1  0.00  0.00 

FREE 

1140 

S  0.00  0.00 

FREE 

1150 

FINISHED 

9760 . 00 


«  ’  «  • 


Example  2--Semicircular  Arch 


84.  The  semicircular  arch  shown  in  Figure  15a  is  replaced  by  an  assem¬ 
blage  of  straight  elements  as  shown  in  Figure  15b.  The  accuracy  of  the  anal¬ 
ysis  using  the  replacement  system  depends  on  the  number  of  nodes  (and  ele¬ 
ments)  used  in  the  model.  The  effect  of  the  number  of  nodes  on  the  solution 
is  shown  in  Figure  15c  for  the  moment  reaction  at  the  left  support  (the  maxi¬ 
mum  bending  moment  for  this  structure). 

85.  Input  data  for  a  41-node  model  were  entered  from  a  predefined  data 
file  as  indicated  below.  An  echoprint  of  input  data  is  shown  on  page  51. 

86.  Results  were  output  selectively  to  the  terminal  as  shown  on 
pages  52-55. 


Input  Data  File  for  Example  2 


1000  EXAMPLE  2  -  SEMICIRCULAR  ARCH 

1010  41  NODE  SOLUTION 

1020  GEOMETRY  F 

1030  1  -15.0  0.0  CURVE  0.0  0.0 

1040  21  0.0  15.0  CURVE  0.0  0.0 

1050  41  15.0  0.0 

1060  PROPERTIES  I  K 

1070  1  41  2000.0  100.0  1000.0  100.0  1000.0 

1080  LOADS  F  F 

1090  D  Y  1  -100.0  21  -100.0 

1100  FIXED  I 

1110  1  0.0  0.0  0.0 

1120  41  0.0  0.0  0.0 

1130  FINISHED 
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Program  Execution  for  Example  2 

PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
BATE:  02/1S/83  TIME:  03(48:35 

ARE  INPUT  DATA  TO  BE  READ  FROM  TERMINAL  OR  FILE7 
ENTER  TERMINAL'  OR  'FILE 

I!  F 

ENTER  INPUT  FILE  NAME  (4  CHARACTERS  MAXIMUM) 
i:  CBNT2 

INPUT  COMPLE1E 

BO  YOU  WANT  INPUT  DATA  ECHOPRINTEB  TO  YOUR 
TERMINAL*  TO  A  FILE*  TO  BOTH*  OR  NEITHER? 

ENTER  'TERMINAL'*  ’FILE'*  'BOTH'*  OR  'NEITHER' 

I  T 


PROGRAM  CBNT  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORT- 
DATE!  02/15/83  TIME!  08:48:5? 

I. --INPUT  DATA 


1 .--HEADING 


EXAMPLE  2  -  SEMICIRCULAR  ARCH 
41  NODE  SOLUTION 


2. --GEOMETRY  DATA 


NODE 

NODE  COORDI 

NATES 

SEGMENT 

CENTER 

COORDINATES 

NO. 

X 

Y 

TYPE 

X 

Y 

(FT) 

(FI  ) 

(FT  > 

'FT  ) 

1 

-15.00 

0.00 

CURVE 

0.00 

0 . 00 

21 

0.00 

15.00 

CURVE 

0,00 

0.00 

41 

15.00 

0.00 

3.-- 

SECTION 

PROPERTIES 

-"-atuiUM  r  KurtK  i  i  tb-- 

———————  •• 

NODE 

NUMBER 

MODULUS  OF 

< - ST 

ART- . - 

< - stop  — . 

START 

STOP 

ELASTICITY 

AREA 

INERTIA 

AREA 

INERTIA 

(K3I ) 

( SOI N ) 

(IN*»4> 

(SO  IN) 

•:iN**4) 

1 

41 

2.00E+03 

100.00 

1000.00 

100.00 

1000.00 

4. --LOAD  DATA 

4. A. --CONCENTRATED  LOADS 
NONE 


4. B. --DISTRIBUTED  LOADS 

LOAD  < - START - >  < - STOP - IN¬ 
DIRECT  NODE  LOAD  NODE  LOAD 

NO.  • P / F  T }  NO.  (P/FT) 

Y  1  -100.00  23  -100.00 


5. --FIXED  SUPPORT  DATA 

SUPPORT  SPECIFIED  DISPLACEMENTS 


NODE 

X-DISP. 

Y-niSP. 

ROTATION 

NO. 

(IN) 

(IN) 

(RAD) 

1 

0.00 

0,00 

0.00 

41 

0.00 

0.00 

0.00 

4. --LINEAR  SPRING  DATA 
NONE 


7. --NONLINEAR  SPRING  DATA 
NONE 


I  >  N 


nq  YOU  WANT  TO  EDIT  INPUT  DATA?  ENTER  'YES'  OR  'NO' 

h 

DO  YOU  WANT  TO  CONTINUE  SOLUTION?  ENTER  YES'  OR  ' NO ' 

I  >  Y 

DO  YOU  WANT  A'lAl.  FORCE  EFFECTS  ON  SENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>N 


SOLUTION  COMPLETE 

DO  YOU  WANT  RESULTS  WRITTEN  TO  YOUR  TERMINAL  .  TO  A  FILE.  OR  BOTH? 
ENTER  'TERMINAL'.  'FILE'.  OR  'BOTH 

I:  T 

FOLLOWING  OUTPUT  DATA  ARE  AVAILABLE ! 

SECTION 


NUMBER  SECTION  CONTENTS 

1 . NODE  COORDINATES  AND  DISPLACEMENTS 

2. . ELEMENT  END  FORCES 

3 . REACTIONS  AT  FIXED  SUPPORTS 


COMPLETE  OUTPUT  CONTAINS  APPROXIMATELY  ?9  LINES. 
OUTPUT  CONTROL 

ENTER  'ALL'.  'SELECTIVE'.  'FINISHED'.  OR  HELP'. 

IIS 

ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT. 

i  l 

ENTER  DESIRED  OUTPUT  UNITS. 

'INCHES'  OR  'FEET'  AND  'FOUNDS'  OR  'KIPS'.  OR  DEFAULT' 

ID 


(Continued) 


Program  Execution  for  Example  2  (Continued) 


INDICATE  AMOUNT  OF  OUTPUT  DESIRED. 

ENTER  All.'.  RANGE  OF  NODES  RETUEEN  1  AND  41.  OR  'HELP'. 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  REAM-COLUMNS  WITH  NONLINEAR  Ml  PORTS 
date;  02/15/83  TIMES  08:55:22 

HEADING 

EXAMPLE  2  -  SEMICIRCULAR  ARCH 
41  NODE  SOLUTION 


NODE  COORDINAT  F  8  AND  DISPLACEMENTS  IN  GLOBAL  COORDINATE  SiSIEM 


MODE 

NODE  100RD1 NAILS 

- node 

D1SPLACEMF.N  13-  - - 

NO  . 

X  (IN) 

Y  (IN) 

(IN) 

Y  (IN' 

ROT  (RAD) 

1 

- 180.00 

0.00 

0 » 

0  , 

0 » 

2 

-179.45 

14.12 

2.2781  -0 

-1 . 5  9  7  E -04 

-3. 109E-04 

3 

-17/ .78 

28.16 

8.39SE-03 

-9 . 506E-04 

5 ; 307E-04 

4 

-175.03 

42.02 

1 . 729E-02 

-2 . 784E-03 

7  2  3  9  E  -  0  4 

J 

-171.19 

.62 

7 . 795E-02 

-5. 85 1C -03 

-8.353E-0-1 

6 

-166.30 

68 . 88 

7.943E-02 

-1 .015E-0? 

-8.901E-04 

*7 

-160.33 

31.72 

5 . 091E-02 

•1 .5491 -02 

- -8.938E-04 

-153.48 

94.05 

6. 169E  -02 

-2. 159E  -02 

-8. 523E  ')4 

? 

-145 : 62 

105.30 

7 . 1 24E--02 

- 2.B03E-02 

-7.71SE-04 

10 

-136.87 

116.90 

7.918E-02 

-3 . 435E-02 

-6 . 586E-04 

li 

-127.28 

127.28 

3 . 528E-02 

-4 . 005E-02 

-5. l«3E-04 

12 

-116.90 

136.67 

8. 943E-02 

-4.466E-0? 

-3.607E- a 4 

13 

-  105.80 

145.62 

9. 186E-02 

-  A . 775E-02 

-1.896E -04 

14 

-94.05 

153.48 

9 , 263E-02 

-4 . 89T'E-02 

-1 .295E  05 

1  5 

-81.72 

160.38 

9.207E-02 

-4 . 006E-02 

1 .623E-0-) 

16 

-68.88 

166.30 

9.057E-02 

-4 • 491E-02 

3 . 291E-04 

1? 

-55.62 

171.19 

8 . 854E-02 

-3 . 954E-02 

4 . S05E-04 

18 

-42.02 

175.03 

8 . 640E-02 

-3.21  IE-02 

6.099E-0J 

19 

-28.16 

177.78 

8 . 452E-02 

-2 . 293E-02 

7. 104E-04 

20 

-14.12 

179.45 

8.323E-02 

-1 . 246E-02 

7 . 756E-04 

21 

0.00 

180.00 

8 . 274E-02 

-1 . 292E- 03 

7 . 99 1 E -04 

22 

14.12 

179.45 

8.313E-02 

9 . 892E-03 

7 . 770E-04 

23 

28.16 

177.78 

8 . 432E-02 

2 . 040E-02 

7 . 1 32E-04 

24 

42.02 

175.03 

8 . 6 10E-02 

2 . 965E-0? 

6 . 1 4  1  E  -  0  4 

25 

55.6? 

171.19 

8 . 316E-02 

3 . 718E-02 

4 . 360E-04 

26 

68.83 

166.30 

9.012E-02 

4 . 26  7E -02 

3 . 357E-04 

27 

81.72 

160.38 

9. 156E-02 

4.595E-02 

1 .700E-04 

28 

94.05 

153.48 

9.208E-02 

4 . 701E-02 

-4.266F.  -06 

2V 

105.80 

145.62 

9. 131E-02 

4.596E-02 

-1 . 8 0 1 E  04 

30 

116.90 

136.87 

8 . 893E-02 

4 . 304E-02 

- ) . SC7E-04 

31 

127.28 

127.23 

8.475E-02 

3 . 861E-02 

-5.0S9E-04 

32 

136.87 

116.90 

7.868E-02 

3 . 309E-02 

-6. 480E-04 

33 

145.62 

105.80 

7 . 080E-02 

2 . 695E-02 

-7 . 613E-04 

34 

153.48 

94,05 

6. 131E-02 

2 . 068E-02 

-9. 421E-04 

35 

160.38 

81.72 

3 . 060E-02 

1 .474E-02 

-8.842E- 04 

36 

166.30 

68.88 

3.919E-02 

9 . 542E-03 

-8.813E  -04 

37 

171 . 19 

55.62 

2 . 778E-02 

S.786E-  03 

-8 . 277E-04 

38 

175.03 

42.02 

l .719E-02 

2 . 448E-03 

-7. 178C-04 

39 

177.78 

28.16 

8 . 349E-03 

7 . 333E-04 

-5 . 463E-04 

40 

179.45 

14.12 

2 . 266E-03 

5 . 320E-05 

-3 . 085E-04 

41 

130.00 

0.00 

0. 

0. 

0. 

ENTER 

SECTION  NUMBER 

FOR  DESIRED 

OUTPUT. 

(Continued) 


Program  Execution  for  Example  2  (Continued) 


K*V 


ENTER  DESIRED  OUTPUT  UNITS. 

'INCHES'  OR  'FEET'  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 


I  >  D 


INDICATE  AMOUNT  OF  OUTPUT  DESIRED. 

ENTER  'ALL'.  RANGE  OF  ELEMENTS  BETWEEN  1  AND 
I>1  TO  5 


40.  OR  'HELP' 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  Bl  AM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE5  02/15/83  TIME:  08:56:14 


HEADING 

EXAMPLE  2  -  SEMICIRCULAR  ARCH 
41  NODE  SOLUTION 


ELEMENT  END  FORCES  IN  ELEMENT  COORDINATE  SYSTEM 


EL. 

NO. 


X-FORCE 
<  P ) 

9 . 922E+02 
9 . 364E+02 
8 . 839E+02 
8 . 353E+02 


-LEFT  END--- 
Y-FORCE 
(P) 

7.911E+02 
7 . 491E+02 
7 . 029E+02 
6 . 526E+02 


MOMENT 
(P-IN) 
4 . 930E+04 
3 .  895E+04 
2 . 920E+04 
2 . 010E+04 


X-FORCE 

<P> 

-9.922E+02 
-9 . 364E+02 
-8 . B39E+02 
-8 . 353E+02 


7 . 908E+02  5 . 987E+02  1.171E+04  -7.908E+02 


-RIGHT  END-- 
Y-FORCE 

<P> 

- 6 .  733E  +  02 
-6.313F+02 
-5 . 351E+02 
-5 . 348E+02 
-4 . 810E+02 


MOMENT 
(P-IN) 
.89SE+04 
i  920E+04 
0 1 0E+04 
1 7  IE  +  04 
076E+03 


DO  YOU  WANT  TO  OUTPUT  ANOTHER  RANGE  FOR  THIS  SECTION? 
ENTER  'YES'  OR  'NO'. 


I>Y 


INDICATE  AMOUNT  OF  OUTPUT  DESIRED. 
ENTER  ALL'.  RANGE  OF  ELEMENTS  BETWEEN 
I>19  23 

ERROR  ON  LAST  ENTRY--REENTER . 

I > 1 9  TO  23 


1  AND  40.  OR  'HELP'. 


•r  v 
r.\\j 

19 

6 . 964E+02 

-3  >  269E+02 

-1 . 181E+04 

-6 . 964E  +  02 

4 . 447E+02 

6  •  358E  +  03 

L%V* 

,  20 

7 . 291E+02 

-3 . 887E+02 

-6 . 358E+03 

-7 . 291E+02 

5 . 065E+02 

3 . 239E+01 

Cv\ 

i  21 

7.666E+02 

-4 . 477E+02 

-3 . 239E+01 

-7.666E+02 

4 . 477E+02 

-6 . 295E+03 

i  22 

7.994E+02 

-3.862E+02 

6.295E+03 

-7 . 994E+02 

3 . 862E  +  02 

-1 . 1 75E  +  04 

S3® 

!  23 

8 . 272E+02 

-3 . 223E+02 

1 . 175E+04 

-8.272E+02 

3 . 223E+02 

-1 . 631E+04 

DO  YOU  WANT  TO  OUTPUT  ANOTHER  RANGE  FOR  THIS  SECTION? 
ENTER  'YES'  OR  'NO', 


I>Y 

INDICATE  AMOUNT 

OF  OUTPUT 

DESIRED. 

n 

ENTER  'ALL'.  RANGE  OF  ELEMENTS  BETWEEN 

I >36  TO  41 

ERROR  ON  LAST  ENTRY--RE ENTER. 

I >36  TO  40 

1  AND  40. 

OR  ' HELP '  . 

36 

7 . 070E+02 

5.3691 +02 

-3 . 796E+03 

-7 . 070E+02 

-5 . 369E+02 

1 . 138E+04 

37 

6 . 627E+02 

5 . 907E+02 

-1 .138E  +  04 

-6 . 627E+02 

-5 . 907E+02 

1 .973E  +  04 

38 

6. 143E+02 

6 . 409E+02 

-1 . 973E  +  04 

-6. 143E+02 

-6 . 409E+02 

2 . 879E  +  04 

39 

5.621E+02 

6 . 871E+02 

-2 . 979E  +  04 

-5 . 621E+02 

-6.871E+02 

3 . 850E  +  04 

40 

5.065E+02 

7.291E+02 

-3 . 850E+04 

-5 . 065E+02 

-7 . 291E+02 

4 . 8B1E  +  04 

DO  YOU  WANT  TO  OUTPUT  ANOTHER  RANGE  FOR  THIS  SECTION? 
ENTER  'YES'  OR  'NO'. 


(Continued) 
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Program  Execution  for  Example  2  (Concluded) 


ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT. 

ENTER  DESIRED  OUTPUT  UNITS. 

'INCHES'  OR  'FEET'  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT- 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE:  02/15/BA  TIME:  08!57:38 

HEADING 

EXAMPLE  2  -  SEMICIRCULAR  ARCH 
41  NODE  SOLUTION 


REACTIONS  AT  FIXED  SUPPORTS  IN  GLOBAL  COORDINATE  SYSTEM 
NODE  X-REACT  Y  REACT  MOM-REACT 

NO.  (P)  (P)  ( P - 1 N ) 

1  -7 .  S16E  +  02  1 . 023E+03  4.930E404 

41  -7 . 484E  +  02  4.770E+02  4.881E404 

ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT. 

DO  YOU  UANT  OUTPUT  WITH  DIFFERENT  UNITS'?’ 

ENTER  'YES'  OR  'NO', 

OUTFUT  COMPLETE 

DO  YOU  UANT  TO  EDIT  INPUT  DATA  FOR  THt  PROBLEM 
JUST  COMPLETED?  ENTER  YES  OP  'NO 

DO  YOU  UANT  TO  MAKE  ANOTHER  'CBNTBM'  RUN?  ENTER  'YES'  OR  'NO' 

**********  NORMAL  TERMINATION  ********** 


87.  The  pile  bent  shown  in  Figure  16a  was  modeled  for  analysis  as  shown 
in  Figure  16b.  The  resistance  of  the  soil  was  treated  as  distributed  linearly 
elastic  springs  with  uniform  stiffness  for  lateral  support  and  as  concentrated 
linearly  elastic  springs  acting  parallel  to  the  pile  axis  at  the  base.  Be¬ 
cause  an  approximate  stiffness  matrix  is  used  for  elements  with  distributed 
springs,  the  number  of  nodes  (elements)  used  for  the  sections  of  the  piles 
below  ground  will  influence  the  accuracy  of  the  solution. 

88.  Program  prompts  and  user  responses  are  shown  on  pages  59  and  60. 

Data  were  entered  from  the  predefined  file  listed  on  page  58.  All  output  data, 
including  an  echoprint  of  input  data  and  selective  output,  were  directed  to  a 
file.  The  output  file,  listed  following  normal  termination  of  the  run,  is 
shown  on  pages  61-64. 

89.  Solutions  for  this  structure  using  10  and  20  elements  to  represent 
the  embedded  portions  of  the  piles  exhibited  only  small  differences  (less  than 
1  percent)  in  results  from  those  given  on  pages  61-64. 
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SEE 


Input  Data  File  for  Example  3 


1000 

EXAMF'Lt  3  -  FILE 

BENT  WITH  LINEAR  SFRINi 

1010 

GEOMETRY  I 

1020 

1 

0.  0, 

1030 

16 

15.  180. 

1040 

17 

22.  264. 

1050 

18 

30.75  369. 

1060 

28 

270.8  369 

1070 

29 

279.5  264 . 

1080 

30 

286.5  180. 

1090 

45 

301.5  0. 

1100 

PROPERTIES  I  K 

1110 

1 

18  4000.  707. 

39761.  707.  39761. 

1120 

18 

28  4000.  1452. 

255552.  1452.  255552. 

1130 

28 

45  4000.  707.  39761.  707.  39761. 

1140 

LOADS  F  K 

1150 

C 

17  1.8  0.0  0 

.0 

1160 

C 

18  18.3  0.0 

0.0 

1170 

c 

28  18.3  0.0 

0.0 

1180 

c 

30  1.8  0.0  0 

.0 

1190 

D 

Y  18  -36.5  28 

-47. 

1200 

LINEAR  I  K 

1210 

C 

1  -94.76  2000 

.  0.0 

1220 

C 

45  -85.24  2000 

.  0. 

1230 

D 

Y  1  2.0  16 

2.0 

1240 

D 

Y  30  2.0  45 

2.0 

1250 

FINISH 

Program  Execution  for  Example  3 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BE Ai  COLUMNS  WITH  NONLINEAR  SUPPORTS 
date:  02/10/83  time:  o?:i6:i7 

ARE  INPUT  DATA  TO  BE  READ  FROM  TERMINAL  OR  FILE? 

ENTER  TERMINAL'  OR  FILE' 

i:  f 

ENTER  INPUT  FILE  NAME  <  b  CHARACTERS  MAXIMUM) 

I  CBN T 3 

INPUT  COMPLETE 

DO  YOU  UANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL-  TO  A  FILE-  TO  BOTH.  OR  NEITHER? 

ENTER  'TERMINAL'-  'FILE'.  'BOTH',  OR  'NEITHER' 

I  F 

ENTER  OUTPUT  FILE  NAME  (6  CHARACTERS  MAXIMUM) 

I>CBNT30 

DO  YOU  UANT  TO  EDI  I  INPUT  DATA?  ENTER  'YES'  OR  'NO' 

I.  N 

DO  YOU  UANT  10  CONTINUE  SOLUTION?  ENTER  'YES'  OR  NO' 

I  Y 

DO  YOU  UANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>N 

SOLUTION  COMPLETE 

DO  YOU  UAMT  RESULTS  WRITTEN  TO  YOUR  TERMINAL. 

TO  FILE  ' CBNT30 '  .  OR  BOTH? 

ENTER  'TERMINAL'.  'FILE'.  OR  'BOTH' 

I>F 

FOLLOWING  OUTPUT  DATA  ARE  AVAILABLE: 


I>S 

I>1 


i:  d 


SECTION 

NUMBER  SECTION  CONTENTS 

1  . NODE  COORDINATES  AND  DISPLACEMENTS 

2  . ELEMENT  END  FORCES 

4  . FORCES  IN  CONCENTRATED  LINEAR  SPRINGS 

5  . FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 

COMPLETE  OUTPUT  CONTAINS  APPROXIMATELY  155  LINES, 

OUTPUT  CONTROL 

ENTER  'ALL'.  'SELECTIVE'.  'FINISHED'.  OR  'HELP'. 


ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT, 


ENTER  DESIRED  OUTPUT  UNITS. 

'INCHES'  OR  'FEET'  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 


INDICATE  AMOUNT  OF  OUTPUT  DESIRED. 

ENTER  'ALL'.  RANGE  OF  NODES  BETWEEN  1  AND  45.  OR  'HELP'. 


I >17  TO  2? 


DO  YOU  WANT  TO  OUTPUT  ANOTHER  RANGE  FOR  THIS  SECTION? 
ENTER  'YES'  OR  NO'. 

i:>n 

ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT. 

I  >2 

ENTER  DESIRED  OUTPUT  UNITS. 

INCHES'  OR  'FEET'  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

I  >  F  K 


(Continued) 


59 


Program  Execution  for  Example  3  (Concluded) 


INDICATE  AMOUNT  OF  0U1PUT  DESIRED. 

ENTER  ' ALL '  .  RANGE  OF  ELEMENTS  BETWEEN  1  AND  44,  OR  'HELP'. 
IMS  TO  28 


DO  YOU  WANT  TO  OUTPUT  ANOTHER  RANGE  FOR  THIS  SECTION? 

ENTER  'YES'  OR  'NO'. 

I,  H 

ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT. 

14 

ENTER  DESIRED  OUTPUT  UNITS. 

INCHES'  OR  'FEET'  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

IM  K 

ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT. 

r  j 

ENTER  DESIRID  OUTPUT  UNITS. 

INCHES'  OR  'FEET'  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

ID 

INDICATE  AMOUNT  OF  OUTPUT  DESIRED. 

ENTER  'ALL'.  RANGE  OF  NODES  BETWEEN  i  AND  44.  OR  'HELP'. 
I  1  TO  16 


DO  YOU  WANT  TO  OUTPUT  ANOTHER  RANGE  FOR  THIS  SECTION? 

ENTER  'YES'  OR  'NO'. 

I>Y 

INDICATE  AMOUNT  OF  OUTPUT  DESIRED. 

ENTER  ALL  .  RANGE  OF  NODES  BETWEEN  1  AND  44.  OR  HELP'. 
I  SO  TO  44 


DO  YOU  WANT  TO  OUTPUT  ANOTHER  RANGE  FOR  THIS  SECTION? 

ENTER  'YES'  OR  'NO'. 

IN 

ENTER  SECTION  NUMBER  FOR  DESIRED  OUTPUT. 

I  F 

DO  YOU  WANT  OU1PUT  WITH  DIFFERENT  UNITS? 

ENTER  'YFS'  OR  NO', 

T  N 

OUTPUT  COMPLETE 

DO  YOU  WANT  TO  EDIT  INPUT  DAIA  FOR  THE  PROBLEM 
JUST  COMPLETED?  ENTER  'YES'  OR  >10' 

r  n 

DO  YOU  WANT  TO  MAKE  AN01T.R  ' CBNTBM '  RUN?  ENTER  'YES'  OR  'NO' 

i:  n 

mtmu*  NORMAL  TERMINATION  ********** 

.440  CP  SECONDS  EXECUTION  TIME 


Output  File  for  Example  3 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
date:  02/15/83  time:  o?:i6:so 

I. --INPUT  DATA 

1  .--HE,'/  ING 

EXAMPLE  3  -  PILE  BENT  WITH  LINEAR  SPRING  SUPPORTS 


2. --GEOMETRY  DATA 


NODE 

NODE  COORDINATES 

SEGMENT 

CENTER 

COGRDINA 

NO. 

X 

V 

TYPE 

X 

« IN) 

(  IN) 

(IN) 

(I 

1 

0.00 

0.00 

STRAIGHT 

16 

15.00 

180.00 

STRAIGHT 

17 

22 . 00 

264.00 

STRAIGHT 

18 

30.75 

369.00 

STRAIGHT 

28 

270.80 

369.00 

STRAIGHT 

29 

279.50 

264.00 

STRAIGHT 

30 

286.50 

130.00 

STRAIGHT 

45 

301.50 

0.00 

3.-- 

SECTION 

PROPERTIES 

< - 

- SECTION 

PROPERTIES 

- > 

NODE 

NUMBER 

MODULUS  OF 

< - 

START - > 

STOP - >• 

START 

STOP 

ELASTICITY 

AREA 

INERTIA 

AREA 

INERTIA 

(KSI  > 

<  SQIN ) 

(IN**4> 

(SQIN) 

( IN«*4> 

1 

18 

4.00E+03 

707.00 

39761.00 

707.00 

39761.00 

18 

28 

4.00E+03 

1452.00 

255552.00 

1452.00 

255552.00 

28 

45 

4 . OOE+03 

707.00 

39761.00 

707. OC 

39761.00 

4 . --LOAD  DATA 

4. A. —CONCENTRATED  LOADS 


CONCENTRATED 

LOADS 

NODE 

X-LOAD 

Y-LOAD 

COUPLE 

NO. 

<K) 

<  K ) 

(K-FT) 

17 

1.80 

0.00 

0.00 

18 

18.30 

0.00 

0.00 

28 

18.30 

0.00 

0.00 

30 

1.30 

0,00 

0,00 

DISTRIBUTED 

LOAD 

LOADS 

DIRECT 

NODE 

LOAD 

NODE 

LOAD 

NO. 

<  K/FT  > 

NO. 

(K/FT) 

Y 

18 

-36.50 

26 

-47.00 

5. --FIXED  SUPPORT  DATA 
NONE 


(Continued) 
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6. --LINEAR  SPRING  DATA 

6. A. --CONCENTRATED  LINEAR  SPRINGS 


SPRING 

SPRING  STIFFNESSES 

NODE 

ANGLE 

TRANSLATION 

ROTATION 

NO. 

(DEG) 

(K/I  ) 

(K-I/RAD) 

1 

-94.76 

2000.00 

0.00 

45 

-35.24 

2000.00 

0.00 

B. --DISTRIBUTED 

LINEAR 

SPRINGS 

SPRING 

< - 

--START - > 

- 

STOP - > 

DIRECT 

NODE 

STIFFNESS 

NODE 

STIFFNESS 

NO. 

(K/SGIN) 

NO. 

(K/SOIN) 

Y 

1 

2.000E+00 

16 

2 .  OOOE+OO 

Y 

30 

2 .  OOOE+OO 

45 

2. OOOE  +  OO 

NONLINEAR  SPRING 

DATA 

NONE 


(Continued) 
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PROGRAM  CBNTBM  -  ANALYSIS  OF  2B  BEAM-COLUMNS  WITH  NONLINEAR  SUFFORTS 
date:  02/15/83  time:  09 : i 7 : 30 

HEADING 

EXAMPLE  5  -  PILE  BENT  WITH  LINEAR  SPRING  SUPPORTS 


NODE 

COORDINATES  AND 

DISPLACEMENTS  IN  GLOBAL 

COORDINATE 

SYSTEM 

NODE 

NODE  COORDINATES 

- NODE 

DISPLACEMF  NTS - > 

NO. 

X  (IN) 

Y  (IN) 

X  (IN) 

Y  (IN) 

ROT  ( R  h  EO 

17 

22 . 00 

264.00 

2.401E-01 

-2.447E-01 

- 1 . 677E-03 

18 

30.75 

369.00 

4.1 73E-01 

-2 . 735E-01 

-1 . 679E- 03 

19 

54.76 

369.00 

4 . 171E-01 

-3. 130E-01 

-1 .576E-03 

20 

78.76 

369.00 

4 . 169E-01 

-3.478E-01 

-1 . 302E-03 

21 

102.77 

369.00 

4. 167E-01 

-3 . 745E-01 

-9.014E-04 

22 

126.77 

369.00 

4. 165E-01 

-3.904E-01 

-4.173E-04 

23 

150.73 

369.00 

4 . 163E-01 

-3. 942E-01 

1 .046E  04 

24 

174.78 

369.00 

4. 162E-01 

-3 . 855E-01 

6. 179E-04 

25 

198.79 

369.00 

4. 160E-01 

-3.650E-01 

1 .075E-03 

26 

222.79 

369.00 

4.158E-01 

-3.347E-01 

1 . 427E-03 

27 

246.80 

369.00 

4 . 1 56E-01 

-2.977E-01 

1 .624E-03 

28 

270.80 

369.00 

4. 154E-01 

-2.584E-01 

1.61 4F- 03 

29 

279.50 

264.00 

4 . 226E-01 

-2 . 405E-01 

-1 .  143E-03 

PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
date:  02/15/83  time:  09:13:01 

HEADING 

EXAMFLE  3  -  PILE  BENT  UITH  LINEAR  SPRING  SUPPORTS 


ELEMENT  END  FORCES  IN  ELEMENT  COORDINATE  SYSTEM 

< . . --LEFT  END - >  < - RIGHT  END-™ 


EL. 

X-FORCE 

Y-FORCE 

MOMENT 

X-FORCE 

Y-FORCE 

MOMENT 

NO. 

(K) 

(K) 

(K-FT) 

(K> 

(K) 

(K-FT) 

18 

4.61 3E+01 

3 . 762E+02 

-1 .506E+01 

-4 . 613E+01 

-  3 . 022E+02 

6 . 939E+02 

19 

4.613E+01 

3 . 022E+02 

-6 « 939E+02 

-4.613E+01 

-2 . 260E+02 

1 . 223E+03 

20 

4.613E+01 

2 . 260E+02 

-1 .223E+03 

-4 . 613E+0 1 

-1 .477E+02 

1 .597E  +  03 

21 

4 . 61 3E+01 

1 . 477E+02 

-1 .597E+03 

-4.61 3E+01 

-6 . 736E+01 

1 .813E+03 

22 

4 . 613E+01 

6 . 736E+01 

-1 .813E+03 

-4 . 613E+01 

1.511EF01 

1 .865EF03 

23 

4.613E+01 

-1.5UE+01 

-1 .865E+03 

-4.613E+01 

9 . 968E+01 

1 .750E+03 

24 

4 . 613E  +  0 1 

-9.963E+01 

- 1 . 750E+03 

-4.613E+01 

1 .863E+02 

1 .464E  +  03 

25 

4.613E+01 

-1 .863E  +  02 

-1 .464E  +  03 

-4.61 3E+01 

2 . 751E+02 

1 .003E+03 

26 

4 . 613E  +  01 

-2 . 751E+02 

-1 .003E+03 

-4.613E+01 

3 . 660EP02 

3. A25E+02 

27 

4.613E+01 

-3.660E+02 

-3.625E+02 

-4 . 613E+01 

4 . 590E+02 

-4 . 623E+02 

28 

4 . 627E+02 

2 . 631E+01 

4 . 623E+02 

-4 , 627E+02 

-2 . 631E+0 1 

-2.313EI02 

(Continued) 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
date:  02/15/83  time:  ov:ia:i6 

HEADING 


EXAMPLE 

3  - 

PILE 

BENT  WITH  LINEAR  SPRING 

SUPPORTS 

FORCES 

IN 

CONCENTRATED  LINEAR  SPRINGS 

- 

-TRANSLATION  SPRING - 1- 

; - ROTATION  SPR 

NODE 

ANGLE 

DEFORMATION  FORCE 

DEFORMATION 

NO. 

<  DEG ) 

(IN)  (K) 

(RAD) 

1 

- 

94.76 

1.885E-01  -3.771E+02 

-1.01 

5E-03  0. 

45 

“ 

85.24 

2.314E-01  -4.629E+02 

-1.88 

5E-03  0. 

PROGRAM 

CBNTBM 

-  ANALYSIS  OF  2D  BEAM-COLUMNS  WITH  NONLINEAR  SUPPOR 

date:  0 

2/15/83 

time:  09: is:: 

HEADING 

EXAMPLE 

3  - 

PILE 

BENT  WITH  LINEAR  SPRING 

supports 

FORCES 

IN 

DISTRIBUTED  LINEAR  SPRINGS  IN  ELEMENT  COORDINATE  SYSTEM 

EL. 

X-SPRING  FORCES  (P/IN) 

Y-SPRING  FORCES  (P/IN) 

NO. 

LEFT 

END  RIGHT  END 

LEFT  END 

RIGHT  END 

1 

0. 

0. 

-1 .657E402 

-1 .412E102 

2 

0. 

0. 

-- 1 . 412E+02 

-1 . 167c F02 

3 

0. 

0. 

-1 . 167E+02 

-9.219E+01 

4 

0. 

0. 

-9.219E+01 

-6.747EF01 

5 

0. 

0. 

-6 . 747E+01 

--4.246E  +  01 

6 

0. 

0. 

-4 . 246E+01 

- 1 . 705E+01 

7 

0. 

0. 

-1 .705E+01 

8 . 921E+00 

8 

0. 

0. 

8. 92IE+00 

3 . 559E+01 

9 

0. 

0. 

3.559E+01 

6.309E+01 

10 

0. 

0. 

6 . 309E+01 

9. 158EF01 

11 

0. 

0. 

9. 158E+01 

1 .212E+02 

12 

0. 

0. 

1 .212E+02 

1 .519E+02 

13 

0. 

0. 

1 .519E+02 

1 .8  ’9E  +  02 

14 

0. 

0. 

1 .839E+02 

2 . 172E+02 

15 

0. 

0. 

2 . 172E+02 

2.517E+02 

16 

0. 

0. 

0. 

0. 

30 

0. 

0. 

-5. 186E+02 

-4 . 693E+02 

31 

0. 

0. 

•  4 . 693E  +  02 

-4. 194E+02 

32 

0. 

0. 

-4.194E+02 

-3 . 695E+02 

33 

0. 

0. 

-3 . 695E+02 

-3. 198E+02 

34 

0. 

0. 

-3 . 198E+02 

-2 . 705E+02 

35 

0. 

0. 

-2 . 705E+02 

-2 • 218E+02 

36 

0. 

0. 

-2 . 218E  +  02 

-1 . 737E+02 

37 

0. 

0. 

-1 .737E+02 

-1 .263E+02 

38 

0. 

0. 

- 1 . 263E+02 

-7 . 937EF01 

39 

0. 

0. 

-7.937E+01 

-3 . 293E+01 

40 

0. 

0. 

-3.298EF01 

1 .300E+01 

41 

0. 

0. 

1 .300E+01 

5.871E+01 

42 

0. 

0. 

5.871E+01 

1 .042EF02 

43 

0. 

0. 

1 .042E  +  02 

1 • 497EF02 

44 

0. 

0. 

1 . 497E+02 

1 .951E+02 

Example  4 — U-Frame  Channel  L iner 

90.  The  U-frame  channel  liner  and  analytical  model  are  shown  in  Fig¬ 
ures  17a  and  17b.  Because  the  structure  and  loading  are  symmetric,  only  half 
of  the  system  was  analyzed.  The  soil  was  modeled  using  nonlinear  distributed 
springs.  Soil  properties  along  the  vertical  walls  were  assumed  to  vary 
linearly  from  zero  at  the  top  to  the  behavior  shown  in  Figure  17c.  Soil  along 
the  base  was  assumed  to  provide  uniform  resistance  as  shown  in  Figure  17d. 

As  in  all  cases  with  distributed  springs,  the  number  of  nodes  (elements)  will 
influence  the  solution. 

91.  The  input  file  and  echoprint  of  input  data  for  the  model  shown  are 
given  on  pages  67-69. 

92.  Complete  output  for  this  example  is  given  on  pages  70-72. 
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Input  Data  File  for  Example  4 


1000 

EXAMPLE  4  - 

U-FRAME 

1010 

GEOMETRY  F 

1020 

1  0.  20. 

1030 

4  0.0  17. 

1040 

5  0.0  16.5 

1050 

6  0.0  16. 

1060 

21  0.0  l. 

1070 

22  0.0  0.5 

1080 

23  0.0  0.0 

1090 

24  0.5  0.0 

1100 

25  1.0  0.0 

1110 

34  10.  0.0 

1120 

PROPERTIES 

I  K 

1130 

1  34  3830.  144.  172' 

1140 

LOADS  F  P 

1150 

D  Y  5  0.0  22 

-1000. 

1160 

D  Y  24  -1000 

.  34  -H 

1170 

FIXED  I 

1180 

34  0.0  FREE 

0.0 

1190 

NONLINEAR  F 

F 

1200 

D  Y  1  2  0.01 

0.0 

1210 

-10.  1.0 

1220 

0.0  0.0 

1230 

E  23  2  0.01 

100, 

1240 

-10.0  1.0 

1250 

72.0  8.0 

1260 

D  Y  23  2  0.01  1000. 

1270 

c 

o 

o 

© 

r4 

1 

1280 

50.0  0.0 

1290 

E  34  2  001 

1000.0 

1300 

-10.0  0.0 

1310 

50.0  0.0 

1320 

FINISHED 

144.  1728. 
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Echoprint  of  Input  Data  for  Example  4 


PROGRAM  CBNTPM 

DATE!  02/15/83 

\m\ 

>> 

I. --INPUT  DATA 

1 . --  HEADING 

ANALYSIS  OF  20  PE  AN  COLUMNS  WITH  NONLINEAR  SUPPORTS 

time:  •:>?: aim? 


EXAMPLE  A  - 

U-FRAME  CHANNEL 

LINER  WITH 

NONLINEAR  SPRING  SUPPORTS 

2. - -GEOMETRY 

DATA 

NODE 

NODE  COORDINATES 

SEGMENT 

CENTER  COORDINATES 

NO. 

X 

V 

TYPE 

X  Y 

(FT) 

(FT) 

.'FT)  (FT) 

1 

0.00 

20.00 

STRAIGHT 

A 

0.00 

17.00 

STRAIGHT 

5 

0.00 

16.50 

STRAIGHT 

6 

0.00 

16.00 

STKAIGH  1 

21 

0.00 

1.00 

STRAIGHT 

22 

0.00 

.50 

STRAIGHT 

23 

0.00 

0.00 

STRAIGHT 

2  A 

•  50 

0.00 

STRAIGHT 

25 

1.00 

0.00 

STRAIGHT 

3  A 

10.00 

0.00 

3. --SECTION 

PROPERTIES 

- 

- SECTION  FROFERTIES - 

NODE  NUMBER 

MODULUS  OF 

- - START - > 

< . —stop - 

START  STOF 

ELASTICITY 

AREA 

INERTIA 

AREA  INERTIA 

( KS  I  ) 

( SQIN ) 

( IN**A) 

( 30 I N )  (IN**A> 

1  3A 

3.83E+03 

1 A A . 00 

1728.00 

1AA.OO  1728.00 

A  < — LOAD  DATA 


*3 


A. A. --CONCENTRATED  LOADS 
NONE 


A. B. --DISTRIBUTED  LOADS 
LOAD  <- 

DIRECT 


NODE 

NO. 


-START - > 

LOAD 
<  F/FT ) 
0.00 

-1000.00 


< . . ST0F - > 

NODE  LOAD 

NO.  (F/FT) 

22  -1000.00 

3A  -1000.00 


5. --FIXED  SUPFORT  DATA 
SUPPORT 
NODE 
NO. 

3A 


SPECIFIED  DISPLACEMENTS 
X-DISF.  Y-DISP.  ROTATION 

<  IN) 

FRFE 


(IN) 

0.00 


(RAD) 

0.00 


6. --LINEAR  SPRING  DATA 
NONE 


(Continued) 
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Echoprint  of  Input  Data  for  Example  4  (Concluded) 


7  .---NONLINEAR  SPRING  BATA 


7. A. --NONLINEAR  CONCENTRATED  SPRINGS 
NONE 

7. B. --NONLINEAR  DISTRIBUTED  SPRINGS 


SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

NODE 

DISPLACEMENT 

FORCE 

NO. 

(FT) 

(P/FT) 

STARTS  Y 

1 

1 . OOE-02 

0. 

DISP.  COORDS.  -10.00  1.00 
FORCE  COORDS.  0.00  0.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

NODE 

DISPLACEMENT 

FORCE 

NO. 

(FT) 

(P/FT) 

ENDS  Y 

23 

1 .OOE-02 

1 .OOE+02 

DISP.  COORDS.  -10.00  1.00 
FORCE  COORDS.  72.00  8.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

NODE 

DISPLACEMENT 

FORCE 

NO. 

(FT) 

(P/FT) 

STARTS  Y 

23 

1 .00E-02 

1 .00E+03 

DISF.  COORDS.  -10.00  0.00 
FORCE  COORDS.  50.00  0.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

NODE 

DISPLACEMENT 

FORCE 

NO. 

(IT) 

(P/FT) 

ENDS  Y 

34 

1  .OOE-02 

1 .00E+03 

DISP.  COORDS 
FORCE  COORDS 


-10.00  0.00 

50.00  0.00 


PROGRAM  CBNTBM  -  ANALYSIS  OF  2D  BEAM-COLUMNS  UITH  NONLINEAR  SUPPORTS 
DATE:  02/15/83  time:  07142:30 

HEADING 

EXAMPLE  4  -  U-FRAME  CHANNEL  LINER  WITH  NONLINEAR  SPRING  SUPPORTS 


NODE  COORDINATES  AND  DISPLACEMENTS  IN  GLOBAL  COORDINATE  SYSTEM 


NODE 

NODE  COORDINATES 

< - NODE 

DISPLACEMENTS - > 

NO. 

X  <  FT  > 

Y  (FT) 

X  (FT) 

Y  (FT) 

ROT  (RAD) 

1 

0.00 

20.00 

4.676E-02 

7.261E  04 

-2.810E-03 

2 

0.00 

19.00 

4 . 395E-02 

7.261E-04 

-2.810E -03 

3 

0.00 

18.00 

4.114E-02 

7.261E-04 

-2.B09E  -03 

4 

0.00 

17.00 

3 . 833E-02 

7.261E -04 

-2.807E-03 

J 

0.00 

16.50 

3.693E-02 

7 . 261E-04 

-2 . 805E-03 

6 

0.00 

16.00 

3.553E-02 

7.261E-04 

-2.801E-03 

7 

0.00 

15.00 

3.273E-02 

7 . 261E-04 

-2. 78BE-03 

8 

0.00 

14.00 

2 . 995E-02 

7.261E-04 

-2 . 765E-03 

9 

0.00 

13.00 

2.721E-02 

7.261E-04 

-2 . 731E-03 

10 

0.00 

12.00 

2.449E-02 

7.261E-04 

-2.685E-03 

11 

0.00 

11.00 

2. 184E-02 

7.261E-04 

-2.624E-03 

12 

0.00 

10.00 

1 .925E-02 

7.261E-04 

-2 . 548E-03 

13 

0.00 

9.00 

1 .675E-02 

7.261E-04 

-2 . 458E-03 

14 

0.00 

8.00 

1.434E-02 

7.261E-04 

-2.354E-03 

15 

0.00 

7.00 

1 .204E-02 

7 . 261E-04 

-2.236E-03 

16 

0.00 

6.00 

9 . 875E-03 

7.261E-04 

-2. 106E-03 

17 

0.00 

5.00 

7.837E-03 

7 . 261E-04 

-1 .965E-03 

18 

0.00 

4.00 

5.945E-03 

7 . 261E-04 

-1 .816E-03 

19 

0.00 

3.00 

4.208E-03 

7.261E-04 

-1 .65BE-03 

20 

0.00 

2.00 

2.633E-03 

7 • 261E-04 

-1 . 43BE-03 

21 

0.00 

1.00 

1 • 236E-03 

7.261E-04 

-1 .304E-03 

22 

0.00 

.50 

6 . 079F.-04 

7.261E-04 

1 .205E-03 

23 

0.00 

0.00 

3.11 BE-05 

7 . 261E-04 

-1 . 100E-03 

24 

.50 

0.00 

2 . 963E-05 

2.033E-04 

-9.911E-04 

25 

1.00 

0.00 

2 . 807E-05 

-2.650E- 04 

-8 . 823E-04 

26 

2.00 

0.00 

2 . 495E-05 

-1 .042E-03 

-6.747E-04 

27 

3.00 

0.00 

2 , 1 83E-05 

-1 .623E-03 

-4 • 915E-04 

28 

4.00 

0.00 

1 .871E-05 

-2 , 036E-03 

-3.399E-04 

2? 

5 » 00 

0.00 

1 . 559E-05 

-2 . 313E-03 

-2.219E-04 

30 

6.00 

0.00 

1.248E-05 

-2.490E-03 

-1 .356E-04 

31 

7.00 

0.00 

9 . 358E-06 

-2 • 594E-03 

-7 . 64BE-05 

32 

8.00 

0.00 

6.237E-06 

-2 . 650E-03 

-3.878E-05 

33 

9.00 

0.00 

3.118E-06 

-2.677E-03 

-1 .569E-05 

34 

10.00 

0.00 

0. 

-2.684E-03 

0. 

(Continued) 
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Output  for  Example  4  (Continued) 


ELEMENT  END  FORCES  IN  ELEMENT  COORDINATE  SYSTEM 


< - 

--LEFT  END-- 

< - 

-RIGHT  END— ■ 

EL. 

X-FORCE 

Y-FORCE 

MOMENT 

X-FORCE 

Y-FORCE 

MOMENT 

NO  . 

<P> 

(P) 

(P-FT) 

(P) 

(P) 

<  P-FT  ) 

1 

0. 

-1 . 080E-07 

-1 . 155E-07 

0. 

-2 . 000E+01 

6 • 667E+00 

2 

0. 

2 . 000E+01 

-6 . 667E+00 

0. 

-8.000E401 

5. 333E+01 

3 

0. 

8 . 000E+01 

-5 . 333E+01 

0. 

-1 .300E+02 

1 .800E+02 

4 

0. 

1 . 800E+02 

-1 . 300E+02 

0. 

-2 . 450E+02 

2 . 858E  +  02 

0. 

2. 450E+02 

-2 . 858E+02 

0. 

-3. 122E+02 

4 . 253E  +  02 

6 

0. 

3.122E+02 

-4 . 253E+02 

0. 

-4.297E+02 

7 . 982E  +  02 

7 

0. 

4 . 297E+02 

-7 • 982E+02 

0. 

-5 . 247E+02 

1 .278E+03 

8 

0. 

5  <  247E+02 

-1 .278E+03 

0. 

-5.972E+02 

1 .840E  +  03 

9 

0. 

5 .  972E  +  02 

-1 .840E+03 

0. 

-6 . 472E+02 

2.464E+03 

10 

0. 

6.472E+02 

-2 . 464E+03 

0. 

-6 . 747E+02 

3. 127E  +  03 

11 

0. 

6 . 747E+02 

-3 . 127E+03 

0. 

-6 . 797E+02 

3 . 806E+03 

12 

0. 

6 • 797E+02 

-3 • 306E+03 

0. 

-6.622E+02 

4 . 479E  +  03 

13 

0. 

6.622E+02 

-4 . 479E+03 

0. 

-6 . 222E+02 

5. 123E+03 

14 

0. 

6.222E+02 

-5 . 1 23E+03 

0. 

-5 • 597E+02 

5 . 716E+03 

15 

0. 

5.597E+02 

-5.716E+03 

0. 

-4.628E+02 

6.228E+03 

16 

0. 

4.628E+02 

-6 . 228E+03 

0. 

-4 . 045E+02 

6 . 657E  +  03 

17 

0. 

4 . 045E+02 

-6 . 657E+03 

0. 

-4. 158E+02 

7 . 061E+03 

18 

0. 

4 . 158E+02 

-7.061E+03 

0. 

-5 . 007E+02 

7 . 513E  +  03 

19 

0. 

5.007E+02 

-7.513E+03 

0. 

-6.617E+02 

8 . 088E  +  03 

20 

0. 

6.617E+02 

-8 . 088E+03 

0. 

-8 • 995E+02 

8 . 858E+03 

21 

0. 

8.995E+02 

-8 • 858E+03 

0. 

-1 .047E+03 

9 . 350E  +  03 

22 

3 . 725E-09 

1 .047E+03 

-9 . 350E+03 

-3.725E-09 

-1 . 720E+03 

1 .003E+04 

23 

1 .720E+03 

-2 . 980E-08 

-1 .003E+04 

-1 . 720E+03 

2.980E-08 

1 .003E+04 

24 

1 .720E+03 

1 . 118E-08 

-1 .003E+04 

-1 .720E+03 

4.858E+02 

9 . 917E+03 

25 

1 i 720E+03 

-4 . 858E+02 

-9  >  917E  +  03 

-1 .720E+03 

1 . 150E+03 

9 . 067E+03 

26 

1 .720E+03 

-1 . 150E+03 

-9.067E+03 

-1 . 720E+03 

1 .477E+03 

7 » 725E  +  03 

27 

1 .720E  +  03 

-1 .477E+03 

-7 • 725E+03 

-1 . 720E403 

1 .556E+03 

6.192E+03 

28 

1 .720E+03 

-1 .556E+03 

-6 • 192E+03 

-1.720E+03 

1 .464E+03 

4.670E+03 

29 

1.720E+03 

- 1 • 464E+03 

-4 . 670E  +  03 

-  1 .720E+03 

1 .259E+03 

3 . 302E+03 

30 

1 . 720E+03 

-1.259EF03 

-3 • 302E+03 

-1 . 720E+03 

9 . 858F+02 

2.175E+03 

31 

1 .720E+03 

-9 . 858E+02 

-2 . 175E+03 

-1 .720E+03 

6 . 732E+02 

1 .343E+03 

32 

1.720E+03 

-6.732E+02 

-1 .343F.+03 

-1 .720E+03 

3 . 407E+02 

8.350E+02 

33 

1 .720E+03 

-3 • 407E  +  02 

-8.350E+02 

-1.720E+03 

1 .490E-08 

6 . 643E+02 

(Continued) 


Output  for  Example  4  (Concluded) 


n 

*  <  ■ 
S : 


REACTIONS  AT  FIXED  SUPPORTS  IN  GLOBAL  COORDINATE  SYSTEM 
NODE  X-REACT  Y-REACT  MOM-REACT 

NO.  (P)  (P)  <P~FT> 

34  -1.720E+03  1.490E-08  6.643E+02 


FORCES 

IN  DISTRIBUTED 

NONLINEAR  SPRINGS 

IN  ELEMENT  COORDINATE  SYSTEM 

EL. 

X-SPRING  FORCES  (P/FT) 

Y- SPRING  FORCES  (P/FT) 

NO. 

LEFT  END 

RIGHT  END 

LEFT  END 

RIGHT  END 

1 

0. 

0. 

0. 

4 .OOOE+Ol 

*? 

0. 

0. 

4 . OOOE+O 1 

8. OOOE+Ol 

3 

0. 

0. 

8 . OOOE+Ol 

1 .200E+02 

4 

0. 

0. 

1 .200E+02 

1 .400E+02 

5 

0. 

0. 

1 . 400E+02 

1 .600F.  +  02 

6 

0. 

0. 

1 .600E+02 

2. 000E  +  02 

7 

0. 

0. 

2 • OOOE+02 

2 . 400E+02 

8 

0. 

0. 

2.400E+02 

2 . 800E  +  02 

9 

0. 

0. 

2 . 800E+02 

3 . 200E+02 

10 

0. 

0. 

3 . 200E+02 

3 . 600E+02 

11 

0. 

0. 

3 . 600E  +  02 

4. OOOE+02 

12 

0. 

0. 

4. 000E+02 

4 . 400E+02 

13 

0. 

0. 

4 . 400E+02 

4 . 800E+02 

14 

0. 

0. 

4 . 800E+02 

5 . 200E+02 

15 

0. 

0. 

5 . 200E+02 

5 . 652E+02 

16 

0. 

0. 

5 . 652E+02 

6 . 944E+02 

17 

0. 

0. 

6.V4E+02 

8 . 287E+02 

18 

0. 

0. 

8 . 287E+02 

9 . 665E+02 

19 

0. 

0. 

9 . 665E+02 

1 . 106E+03 

20 

0. 

0. 

1 . 106E+03 

1 .244E+03 

21 

0. 

0. 

1 .244E+03 

1 .313E+03 

22 

0. 

0. 

1 . 313E+03 

1 .380E+03 

23 

0. 

0. 

0. 

0. 

24 

0. 

0. 

0. 

1 . 325E+02 

25 

0. 

0. 

1 . 325E+02 

5 . 209E+02 

26 

0. 

0. 

5 . 209E+02 

8. 113E+02 

27 

0. 

0. 

8 . 1 13E  +  02 

1 .018E+03 

28 

0. 

0. 

1 . 018E+03 

1 .  157E+03 

29 

0. 

0. 

1 . 157E+03 

1 . 244E  +  03 

30 

0. 

0. 

1 . 244E  +  03 

1 .297E+03 

31 

0. 

0. 

1 .297E+03 

1.325E+03 

32 

0. 

0. 

1 .325E+03 

1 .338E+03 

33 

0. 

0. 

1 .338E+03 

1 .342E  +  03 

Destroy  this  report  when  no  longer.needed.  Do  not  return 
it  to  the  originator. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated 
by  other  authorized  documents. 


The  contents  of  this  report  are  not  to  be  used  for 
advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of 
such  commercial  products. 
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